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Abstract 

Simulations  are  conducted  with  a  cloud-resolving  numerical  model  to  examine  the 
transformation  of  a  weak  incipient  mid-level  cyclonic  vortex  into  a  tropical  cyclone. 
15  Results  demonstrate  that  two  distinct  pathways  are  possible  and  that  development  along  a 
particular  pathway  is  sensitive  to  model  physics  and  initial  conditions.  One  pathway 
involves  a  steady  increase  of  the  surface  winds  to  tropical  cyclone  strength  as  the  radius 
of  maximum  winds  gradually  decreases.  A  notable  feature  of  this  evolution  is  the  creation 
of  small-scale  cyclonic  vorticity  anomalies  by  deep  convective  towers  and  subsequent 
20  merger  and  convergence  by  the  low-level  secondary  circulation.  The  second  pathway  also 
begins  with  a  strengthening  low-level  circulation,  but  eventually  a  significantly  stronger 
mid-level  circulation  develops.  Cyclogenesis  subsequently  occurs  when  a  small-scale 
surface  concentrated  vortex  abruptly  forms  near  the  center  of  the  larger-scale  circulation. 
The  small-scale  vortex  is  warm  core  throughout  the  troposphere  and  results  in  a  local 
25  surface  pressure  fall  of  a  few  millibars.  It  usually  develops  rapidly  undergoing  a  modest 
growth  to  form  a  small  tropical  cyclone.  Many  of  the  simulated  systems  approach  or 
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50  1  Introduction 


The  problem  of  how  a  tropical  disturbance  transforms  into  a  tropical  cyclone  (TC)  has 
been  an  active  field  of  investigation  for  over  fifty  years.  Numerical  modeling  has  become 
an  increasingly  useful  research  tool,  as  more  powerful  computer  processing  and  the 
development  of  advanced  cloud  models  enable  more  accurate  simulation  of 
55  microphysical  processes,  air-sea  interaction  and  radiative  exchange.  Relatively  recent 
studies  that  have  attempted  to  resolve  clouds,  rather  than  parameterize  their  effects,  show 
significant  ability  to  reproduce  observed  characteristics  of  mature  TCs  (Liu  et  al.  1997; 
Braun  2002;  Rogers  et  al.  2003). 

Hendricks  et  al.  (2004,  hereafter  H04)  extended  the  genesis  study  of  Davis  and  Bosart 
60  (2001)  by  analyzing  a  high-resolution  (3  km  horizontal  grid  spacing)  near-cloud¬ 

resolving  numerical  simulation  of  Hurricane  Diana  (1984)  with  the  fifth  generation 
Pennsylvania  State  University-National  Center  for  Atmospheric  Research  (PSU-NCAR) 
Mesoscale  Model  (MM5).  The  numerical  model  was  initialized  with  an  analysis  of 
observed  data  at  a  stage  when  the  system  that  would  evolve  into  Hurricane  Diana  was  a 
65  subtropical  disturbance.  Results  suggested  that  convective  plumes  with  intense  vertical 
vorticity  in  their  cores  were  the  preferred  convective  structures.  The  term  Vortical  Hot 
Tower  (VHT)  was  coined  to  describe  these  structures.  H04  suggested  that  VHT’s  play  an 
important  role  in  a  two-stage  evolutionary  process:  (1)  preconditioning  of  the  local 
environment  via  diabatic  production  of  multiple  small-scale  lower  tropospheric  cyclonic 
70  potential  vorticity  (PV)  anomalies,  and  (2)  multiple  mergers  and  axisymmeterization  of 
these  low-level  PV  anomalies.  In  addition  to  the  organizational  process  of  the  PV 
anomalies,  the  cyclogenesis  is  enhanced  by  the  aggregate  diabatic  heating  associated  with 
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the  VHT’s,  which  produce  a  net  influx  of  low-level  mean  angular  momentum  throughout 
the  genesis. 

75  Similar  results  to  H04  were  obtained  in  a  study  by  Montgomery  et  al.  (2006,  hereafter 
M06).  Unlike  H04,  an  idealized  model  initialization  was  used  based  on  observations  from 
the  1990s  that  indicated  that  mesoscale  convective  vortices  (MCVs)  forming  in  the 
stratiform  region  near  mid-levels  of  the  troposphere  in  the  trough  region  of  tropical 
disturbances  are  often  precursors  to  TC  genesis  in  the  deep  tropics  (e.g.  Zehr  1992; 
80  Mapes  and  Houze  1995;  Harr  and  Elsberry  1996;  Harr  et  al.  1996;  Bister  and  Emanuel 
1997;  Ritchie  and  Holland  1997;  Raymond  et  al.  1998;  Reasor  et  al.  2005).  They  used  the 
Regional  Atmospheric  Modeling  System  (RAMS)  to  explore  the  formation  of  the  TC 
surface  circulation  through  a  series  of  idealized  simulations  of  convection  initiated  within 
a  preexisting  MCV.  For  the  majority  of  experiments  the  maximum  wind  speed  of  the 
85  initial  vortex  was  6.6  m  s'1  at  a  radius  of  75  km  and  a  height  of  4  km.  The  mean  Atlantic 
hurricane  season  sounding  from  Jordan  (1958)  was  used  to  initialize  the  model 
thermodynamic  structure.  The  center  of  the  vortex  was  moistened  at  low  levels  and  the 
sea  surface  temperature  was  set  to  a  constant  value  of  29  °C  for  most  of  the  experiments. 
The  microphysics  scheme  included  the  ice  phase  and  a  radiation  scheme  was  used  for  all 
90  the  simulations  conducted.  An  important  ramification  of  this  study,  was  that  it  provided 
support  for  the  view  that  a  relatively  weak  mid-level  MCV  in  a  favorable  thermodynamic 
environment,  over  a  warm  ocean,  would  gradually  undergo  a  metamorphosis  into  a  TC  in 
the  absence  of  other  influences,  such  as  large  scale  forcing,  or  merger  of  two  or  more 
MCV’s  (e.g.  Simpson  et  al.  1997;  Sippel  et  al  2006).  Factors  such  as  large-scale  forcing 
95  and  merger  of  MCV’s  may  increase  the  likelihood  of  genesis  in  some  situations. 
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However,  the  M06  study  suggested  the  possibility  that  an  embryonic  vortex  in  an 
environment  devoid  of  hostile  influences,  such  as  strong  vertical  wind  shear  or  dry  air 
intrusion,  may  gradually  undergo  a  transformation  into  a  TC. 

Recent  work  by  Dunkerton  et  al.  (2009),  broadened  the  earlier  investigation  of  M06  to 
100  include  synoptic  and  sub-synoptic-scale  influences  associated  with  easterly  waves  and 
their  embedded  critical  layers  and  cyclonically  recirculating  flow  within.  Based  on  their 
examination  of  55  developing  cases  in  the  Atlantic  and  eastern  Pacific  sectors  during  the 
peak  of  four  consecutive  hurricane  seasons,  Dunkerton  et  al.  (2009)  developed  the  new 
“marsupial  paradigm”  for  tropical  cyclogenesis  that  explicitly  recognizes  the  intrinsic 
105  multi-scale  nature  of  the  problem.  The  critical  layer  of  a  tropical  easterly  wave,  or  “wave 
pouch”,  was  hypothesized  to  be  important  to  TC  formation  in  three  distinct  ways:  (1) 
wave  breaking  or  roll-up  of  cyclonic  vorticity  and  lower-tropospheric  moisture  near  the 
critical  surface  in  the  lower  troposphere  provides  a  favorable  environment  for  the 
aggregation  of  vorticity  seedlings  for  TC  formation;  (2)  the  cat’s  eye  is  a  region  of 
110  approximately  closed  circulation,  where  air  is  repeatedly  moistened  by  deep  moist 
convection  and  protected  to  some  degree  from  dry  air  intrusions;  (3)  the  parent  wave  is 
maintained  and  possibly  enhanced  by  diabatically  amplified  mesoscale  vortices  within 
the  wave.  The  new  cyclogenesis  model  is  supported  in  part  by  observations  of  a 
developing  Pacific  easterly  wave  (Montgomery  et  al.  2010,  Raymond  and  Lopez  2011) 
115  and  high  resolution  cloud-representing  numerical  simulations  in  both  real-case  and 
idealized  configurations  (Zhang  et  al.  2010,  Montgomery  et  al.  2010,  Wang  et  al.  2010, 
Montgomery  et  al.  2011).  Despite  these  encouraging  scientific  results,  questions  remain 
about  the  nature  of  the  convective  organization  process  within  the  wave  pouch  in  the 
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early  stages  of  a  tropical  disturbance.  Some  of  the  issues  are  exemplified  by  the  two 
different  cyclogenesis  pathways  reported  in  M06  and  Nolan  (2007,  hereafter  N07). 
Another  open  issue  concerns  the  thermodynamic  nature  of  tropical  cyclogenesis 
(Raymond  et  al.  2011,  Smith  and  Montgomery  2012,  Montgomery  and  Smith  2012). 

Some  of  the  specific  issues  arising  with  the  two  studies  of  M06  and  N07  can  be 
summarized  as  follows:  Although,  M06  found  a  strengthening  of  the  mid-level 
125  circulation  often  occurred  in  their  simulations  (e.g.  Fig.  4c  of  M06),  this  mid-level 
strengthening  did  not  appear  to  play  a  crucial  role  in  the  genesis  process.  In  contrast,  N07 
found  a  significantly  different  pathway  to  tropical  cyclogenesis  occurred  in  simulations 
conducted  with  the  Weather  Research  and  Forecasting  (WRF)  model,  one  in  which  a 
mid-level  vortex  was  suggested  to  play  an  important  role.  Some  of  the  experiments  in 
130  N07  were  initialized  with  a  mid-level  vortex,  similar  to  M06.  The  maximum  strength  of 

this  initial  vortex  was  stronger  with  a  maximum  wind  speed  of  10  m  s’1  at  a  radius  of  100 
km  and  a  height  of  3.72  km.  Unlike  M06,  the  initial  vortex  was  not  moistened  above  that 
of  the  surrounding  environment  and  radiation  was  not  included.  After  approximately 
three  days  simulation  time,  a  strengthening,  and  contraction  of  the  mid-level  circulation 
135  occurred  with  azimuthally  averaged  winds  reaching  about  12  ms’1  at  a  radius  of  65  km 
and  at  a  height  of  5  km.  Low  level  winds  increased  also  to  approximately  8  m  s'1.  At  this 
time,  a  small  surface  concentrated  vortex  (SSCV)  formed  suddenly  near  the  center  of  the 
large  scale  circulation,  with  a  radius  of  about  7  km,  and  subsequently  the  system  rapidly 
developed.  NO  7  emphasized  that  genesis  did  not  occur  until  the  inner  core  had  achieved 
140  deep  near-saturation  and  the  mid-level  vortex  had  elevated,  contracted,  and  intensified. 
In  the  hours  before  genesis,  the  intensification  of  the  mid-level  vortex  was  shown  to  lead 
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to  a  large  increase  in  the  (axisymmetric)  efficiency  of  the  conversion  of  latent  heat  energy 
to  the  kinetic  energy  of  the  cyclonic  wind  field,  in  accord  with  the  theory  of  Schubert  and 
Hack  (1980)  and  Hack  and  Schubert  (1986).  In  some  respects,  this  description  of  the 
145  evolution  resembles  the  conceptual  model  proposed  by  Bister  and  Emanuel  (1997)  who 
also  postulate  a  crucial  role  for  the  mid-level  vortex  in  TC  genesis. 

The  current  study  presents  results  of  simulations  conducted  with  a  newer  version  of 
RAMS  than  used  by  M06.  Unlike  the  previous  version  of  RAMS,  the  simulations  with 
the  newer  version  show  more  of  a  tendency  to  produce  a  prominent  mid-level  vortex. 
150  Additionally,  the  newer  version  sometimes  shows  the  subsequent  formation  of  a  SSCV 
very  reminiscent  of  the  results  of  N07.  In  this  study  we  will  refer  to  the  pathway 
described  by  M06  as  pathway  one  and  the  pathway  described  by  N07  as  pathway  two. 
Experiments  are  performed  to  examine  the  sensitivity  of  the  pathway  taken  to  model 
physics  and  initial  conditions.  We  summarize  the  distinctive  features  of  these  different 
155  pathways  as  follows: 

Pathway  One:  Pathway  one  proceeds  by  spin  up  of  cyclonic  surface  winds  that  become 
greater  than  those  of  the  initial  weak  mid-level  vortex  and  remain  stronger  or  of 
comparable  strength  as  the  winds  aloft,  until  they  become  large  enough  that  they  reach 
160  tropical  depression  strength  (defined  herein  as  approximately  12  m  s’1).  Development  is 
characterized  by  a  gradual  decrease  of  the  radius  of  maximum  surface  winds.  Vorticity  is 
concentrated  in  numerous  deep  convective  towers  that  collectively  drive  a  system-scale 
inflow.  Vorticity  gradually  builds  in  the  center  as  the  system-scale  inflow  produces 
increasing  cyclonic  vorticity  and  as  small-scale  cyclonic  vorticity  anomalies  are 
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165  converged  at  low  levels  and  undergo  aggregation.  The  foregoing  process  is  usually  a 
stochastic  process  with  no  prominent  event  leading  to  a  sudden  fall  in  the  minimum 
surface  pressure. 

Pathway  Two:  For  pathway  two  a  mid-level  vortex  develops  in  excess  of  10  m  s'1  that  is 
170  considerably  stronger  than  the  typical  MCY  associated  with  the  stratiform  region  of  a 
Mesoscale  Convective  System  (MCS).  Its  size  is  usually  smaller  than  the  vortex  used  to 
initialize  the  model,  with  a  radius  of  maximum  winds  ranging  from  approximately  25  to 
100  km,  for  the  experiments  discussed  herein.  Although  the  surface  winds  may 
strengthen  prior  to  the  development  of  the  strong  mid-level  vortex,  with  a  concomitant 
175  decrease  in  the  radius  of  maximum  tangential  winds  similar  to  pathway  one,  the 
intensification  of  the  surface  winds  ends  once  the  strong  mid-level  vortex  is  established. 
The  system  development  is  then  followed  by  the  sudden  formation  of  a  SSCY  near  the 
center  of  the  larger-scale  circulation.  The  formation  of  the  SSCV  is  accompanied  by  a 
local  decrease  of  surface  pressure  of  a  few  millibars.  Once  formed,  the  SSCV  becomes 
180  the  focus  of  a  strengthening  cyclonic  circulation  that  grows  in  size,  and  often  develops 
rapidly  into  a  small  TC. 

Although  the  newer  version  of  RAMS  behaves  differently  than  the  older  version  in 
that  it  is  more  likely  to  show  development  along  pathway  two,  we  should  be  careful  not 
185  to  draw  definitive  conclusions  at  this  stage  because  results  demonstrate  a  significant 
sensitivity  to  model  physics  and  initial  conditions.  Further  refinements  to  model 
parameterizations  could  lead  to  a  different  picture  of  the  likelihood  of  a  particular 
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pathway  being  favored.  The  solution  to  this  initial  value  problem  appears  to  have 
similarities  to  a  bifurcation  in  that  small  changes  to  parameters  seem  to  lead  to  a  sudden 
190  qualitative  change  in  the  behavior.  The  strength  of  the  mid-level  circulation  just  prior  to 
genesis  for  systems  developing  along  pathway  two  in  the  simulations  of  N07  and  those 
discussed  in  this  paper  are  notably  large,  often  in  excess  of  12  m  s'1.  Observations  of 
MCVs  that  would  support  the  existence  of  such  strong  wind  speeds  prior  to  genesis  are 
limited.  Furthermore,  we  are  unaware  of  any  observed  cases  that  clearly  show  the  abrupt 
195  formation  of  a  S SC V. 

As  far  as  we  are  aware,  the  rather  remarkable  transformations  that  take  place  along 
pathway  two  have  been  only  demonstrated  so  far  by  N07  and  the  results  presented  in  this 
modeling  study.  Other  numerical  modeling  studies  by  H04,  M06,  Wang  et  al.  (2010)  and 
Braun  et  al.  (2010),  do  not  show  development  along  pathway  two.  The  initial  conditions 
200  used  in  the  studies  by  N07  and  herein,  are  idealized;  Large  scale  forcing,  complex 
interactions  between  multiple  mesoscale  vortices,  vertical  wind  shear,  and  development 
within  a  tropical  wave  could  all  have  significant  influences  on  the  pathway  taken  to 
tropical  cyclogenesis.  Nevertheless,  it  is  interesting  that  there  are  two  distinct  pathways  to 
tropical  cyclogensis  that  are  found  in  numerical  modeling  simulations  even  with  simple 
205  idealized  initial  conditions,  and  it  is  important  to  understand  why  this  occurs. 

As  discussed  by  Montgomery  and  Smith  (2010),  the  results  of  N07,  raises  several 
questions:  What  processes  lead  to  the  intensification  and  contraction  of  the  mid-level 
vortex?  Although  saturation  of  the  inner  core  and  strengthening  of  the  mid-level  vortex 
occur  prior  to  the  formations  of  the  SSCV,  are  there  other  factors  that  influence  its 
210  development?  The  small  scale  vortex  subsequently  intensifies  and  grows  into  the  primary 
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circulation,  but  how  this  happens  requires  explanation.  N07  considers  the  formation  of 
the  SSCV  to  be  the  time  of  tropical  cyclogenesis.  Yet,  given  its  small-scale  feature,  and 
given  the  weak  larger-scale  surface  circulation,  the  classification  of  the  system  as  a  TC  at 
this  stage  seems  questionable  to  us.  Moreover,  the  metamorphosis  process  described  by 
215  N07,  which  is  in  sharp  contrast  to  the  one  found  in  the  earlier  studies  by  H04  and  M06, 

may  not  be  the  typical  pathway  along  which  most  TCs  form.  For  our  simulations 
discussed  in  this  paper  that  evolve  along  pathway  two,  the  simulated  TC  usually  remains 
very  small,  which  suggests  that  it  is  unlikely  to  be  the  way  that  most  TCs  form. 
Nevertheless,  it  may  turn  out  to  be  an  important  mechanism  for  the  formation  of  small 
220  TCs,  which  are  hard  to  forecast.  Also,  it  is  important  for  hurricane  modelers  and 
forecasters  to  be  aware  that  such  a  pathway  may  exist,  so  that  they  recognize  it  if  it 
occurs  in  their  simulations  or  observations.  This  exploratory  study  seeks  to  further  our 
understanding  of  why  a  particular  pathway  may  be  favored  and  elucidate  some  of  the 
factors  responsible  for  evolution  along  pathway  two. 

225  An  outline  of  the  remaining  paper  is  as  follows:  In  section  2,  we  discuss  the  numerical 
model,  and  the  initial  conditions  used  for  the  experiments.  In  section  3,  results  are 
presented.  Simple  statistics  are  compiled  for  the  experiments  so  that  their  general 
behavior  can  be  compared  and  contrasted.  Individual  cases  are  then  analyzed  in  more 
detail.  This  is  followed  by  an  examination  of  the  transformations  that  take  place  during 
230  pathway  two,  in  particular  the  processes  that  appear  to  be  important  in  the  development 
of  a  prominent  mid-level  vortex  and  subsequently  the  abrupt  formation  of  a  SSCV. 
Results  are  discussed  in  section  4  and  conclusions  are  presented  in  section  5. 
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2  Model  description  and  initial  conditions 

235 

2.1  Cloud  model 

The  model  is  a  more  recent  version  of  the  one  used  by  M06.  It  is  very  similar  in  construct 
although  numerous  small  changes  were  made  to  the  code.  The  main  features  are  that  it  is 
240  a  nonhydrostatic  numerical  modeling  system  comprising  time-dependent  equations  for 
velocity,  non-dimensional  pressure  perturbation,  ice-liquid  water  potential  temperature 
(Tripoli  and  Cotton  1981),  total  water  mixing  ratio  and  cloud  microphysics.  The 
microphysics  scheme  has  categories  for  cloud  droplets,  rain,  pristine  ice  crystals,  snow, 
aggregates  and  hail  (Walko  1995).  The  surface  parameterization  of  heat,  vapor  and 
245  momentum  fluxes  is  based  on  the  Louis  (1979)  scheme.  The  roughness  length  over  water 
is  dependent  on  the  surface  wind  speed  according  to  the  relation  derived  by  Charnock 
(1955).  It  includes  the  longwave  and  shortwave  scheme  developed  by  Harrington  (1997), 
which  is  employed  for  most  of  the  experiments.  A  simpler  scheme  developed  by  Chen 
and  Cotton  (1987)  is  an  option  in  RAMS  and  is  used  for  a  sensitivity  experiment.  Both 
250  schemes  allow  for  interactions  of  the  radiation  with  cloud  hydrometeors.  A  standard  first- 
order  sub-grid  scale  turbulence  scheme  developed  by  Smagorinsky  (1963)  is  used  with 
modifications  by  Lilly  (1962)  and  Hill  (1974)  that  enhance  diffusion  in  unstable 
conditions  and  reduces  diffusion  in  stable  conditions.  RAMS  utilizes  the  two-way 
interactive  multiple  nested  grid  scheme  developed  by  Clark  and  Farley  (1984).  The 
255  radiative  boundary  condition  described  by  Klemp  and  Wilhelmson  (1978)  is  used  at  the 
lateral  boundary  of  the  coarse  grid.  A  Rayleigh  friction  layer  is  included  at  upper  levels. 
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The  values  of  the  model  parameters  employed  remain  the  same  as  in  the  previous  study 
(see  the  appendix  of  M06),  except  that  the  value  of  Charnocks  constant  is  reduced  from 
0.018  to  0.016. 

260  For  all  the  numerical  experiments  three  grids  are  used.  The  horizontal  grid  increments 
are  24,  6  and  2  km,  with  (x,y,z)  dimensions  of  101x101x33,  110x1110x33,  161x161x33, 
respectively  for  the  small  vortex  simulations  and  101x101x33,  1 10x1 10x33,  201x201x33, 
respectively  for  the  large  vortex  simulations.  Each  grid  is  centered  within  the  next 
coarsest  grid.  The  vertical  grid  increment  is  200  m  and  gradually  stretched  with  height  to 
265  the  top  of  the  domain  at  z=23.4  km.  The  depth  of  the  Rayleigh  friction  layer  is  6  km. 

2.2  Initial  environmental  temperature  and  moisture  profiles 

Similarly  to  M06,  the  temperature  structure  is  the  mean  Atlantic  hurricane  season 
270  sounding  of  Jordan  (1958).  The  initial  moisture  profile  is  the  same  as  used  for  the 
reduced  CAPE  experiments  in  M06,  referred  to  as  B2  and  B3  in  that  study.  This  profile 
has  low  level  moisture  reduced  from  the  Jordan  sounding  by  a  maximum  of  2  g  kg'1  at 
the  surface.  The  reduced  moisture  of  this  profile  may  be  more  representative  of  the 
environment  surrounding  a  TC,  which  would  be  expected  to  have  a  moister  core.  It  also 
275  has  the  advantage  of  producing  a  more  focused  development,  so  that  there  is  less 
convective  activity  in  the  coarser  grids  that  only  poorly  resolve  convective  scales. 

2.3  Experiments 
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The  details  of  the  procedure  used  to  initialize  the  MCV,  is  discussed  by  M06. 
Maximum  tangential  winds  for  all  the  experiments  are  at  a  height  of  4  km  above  sea 
level.  Experiments  are  broadly  grouped  into  the  following  categories:  (1)  Radiation 
versus  no  radiation  Simulations  are  conducted  either  with  the  Harrington  scheme 
activated  or  with  it  deactivated.  (2)  Moist  versus  dry.  For  most  experiments  the  initial 
285  vortex  is  moistened  below  8  km  and  for  a  radii  less  than  the  radius  of  maximum  winds 
(RMW),  to  85%  of  saturation  with  respect  to  liquid.  This  moisture  anomaly  is  linearly 
reduced  to  environmental  values  from  the  RMW  to  a  radius  of  25  km  beyond  the  RMW. 
For  the  dry  simulations  there  is  no  vortex  moistening  and  the  moisture  is  horizontally 
homogeneous  at  environmental  values.  (3)  Small  versus  large  vortices.  Experiments  are 
290  conducted  for  two  vortex  sizes;  One  having  a  RMW  of  75  km,  and  one  with  a  RMW  of 
125  km.  (4)  Weak  versus  strong.  Most  of  the  experiments  are  initialized  with  a  relatively 
weak  vortex  having  a  maximum  wind  speed  of  8  m  s'1  at  a  height  of  4  km  above  the 
surface,  and  a  maximum  surface  wind  speed  of  4  m  s’1.  A  few  experiments  were  carried 
out  with  a  stronger  vortex  having  a  maximum  wind  speed  of  12  m  s'1  at  a  height  of  4  km 
295  above  the  surface,  and  a  maximum  surface  wind  speed  of  6  m  s'1.  (5)  Sea  surface 

temperature  (SST).  The  majority  of  experiments  used  an  SST  of  29  °C.  Some  simulations 
were  performed  with  SSTs  of  28  °C  and  27  °C. 

As  an  example  of  an  initial  vortex  Figure  1  shows  the  tangential  wind  speed,  potential 
temperature  perturbation  and  water  vapor  mixing  ratio  for  a  small  weak  vortex  that  has  a 
300  moisture  anomaly.  The  maximum  wind  speeds  are  at  a  height  of  4  km  and  a  radius  of  75 

km.  The  potential  temperature  anomaly  is  a  maximum  at  z=6  km  and  a  minimum  at  z=2 
km  at  the  center  of  the  vortex.  Moistening  to  85%  of  saturation  below  8  km  within  the 
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vortex  increases  the  water  vapor  mixing  ratio  by  approximately  2  g  kg’1  relative  to  the 
environment.  The  effects  of  the  moisture  anomaly  are  included  in  the  buoyancy  and 
305  virtual  potential  temperature  used  to  specify  the  initial  vortex  (M06),  so  that  it  remains 
well  balanced.  These  changes  are  reflected  in  the  potential  temperature  perturbation  field 
particularly  evident  where  the  horizontal  gradient  of  water  vapor  is  strong  between  75- 
100  km.  The  saturated  vapor  pressure  is  a  function  of  temperature  that  is  perturbed  in  the 
vortex  resulting  in  a  water  vapor  mixing  ratio  field  that  is  not  completely  flat  in  regions 
310  that  are  85%  of  saturation. 

Table  1  shows  15  experiments  numbered  2  to  16  conducted  using  the  grouping 
discussed  above.  All  these  experiments  have  ice  microphysics.  Experiment  1  that  is  not 
included  in  this  table  is  the  only  simulation  that  does  not  have  ice  microphysics.  This 
simpler  microphysics  case  is  discussed  first.  It  is  identical  to  Experiment  2  except  it  only 
315  has  cloud  water  and  rain  microphysical  categories.  Experiments  2  to  7  do  not  have 
radiation.  These  simpler  physics  cases  are  discussed  prior  to  experiments  8  to  16  that 
include  radiation.  Additionally,  one  more  simulation,  Experiment  17,  was  conducted 
using  the  Chen-Cotton  radiations  scheme  instead  of  the  Harrington  scheme. 

Experiment  2  is  similar  to  the  reduced  CAPE  simulations  B2  and  B3  of  M06  except 
320  that  the  form  of  the  moist  anomaly  is  different  and  radiation  is  not  included.  Experiment 
3  is  identical  to  Experiment  2  except  the  initial  vortex  is  dry.  Experiment  4  is  similar  to  2 
except  the  initial  vortex  is  large.  Experiments  5  and  6  have  a  lower  SST  of  28  °C  for  a 
small  and  large  vortex,  respectively.  Experiment  7  is  identical  to  experiment  2  except  the 
SST  is  reduced  from  29  to  27  °C.  Experiment  8  is  identical  to  Experiment  2  except  that 
325  radiation  is  activated.  Experiment  9  is  identical  to  8  except  the  initial  vortex  is  dry. 
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Experiments  10  and  11  are  stronger  vortex  cases  for  moist  and  dry  small  vortices 
respectively.  Experiment  12  is  identical  to  8  except  the  vortex  is  large.  Experiment  13  is 
identical  to  experiment  12  except  the  vortex  is  strong.  Experiments  14-16  examine  the 
effects  of  reducing  SST  when  radiation  is  activated.  Finally,  Experiment  17  is  identical  to 
330  Experiment  8  except  that  the  Chen-Cotton  radiation  scheme  is  employed.  Clearly,  many 
more  combinations  are  possible.  The  experiments  that  have  been  chosen  should  give  a 
fairly  comprehensive  view  of  the  affects  of  changing  some  of  the  major  parameters  on  the 
pathway  taken  to  tropical  cyclogenesis. 

335  3  Experimental  results 

3.1  Discussion  of  the  general  behavior 

Table  2  shows  some  general  statistics  for  the  experiments  that  include  the  time  at 
340  which  the  maximum  averaged  tangential  wind  at  the  grid  point  adjacent  to  the  surface 
(z=98m)  reaches  12  m  s'1,  and  the  RMW  at  this  time.  For  systems  developing  along 
pathway  one,  this  wind  speed  is  a  reasonable  criteria  for  concluding  that  tropical 
cyclogenesis  has  occurred.  For  systems  developing  along  pathway  two  that  form  an 
SSCV  this  wind  speed  is  probably  not  a  good  indicator  of  genesis  since  the  larger  scale 
345  circulation  at  the  surface  is  still  weak.  Also  shown  is  the  pathway  taken  to  genesis. 
Designation  as  pathway  two  requires  the  development  of  a  prominent  mid-level  vortex 
and  the  subsequent  abrupt  formation  of  a  resilient  SSCV  that  becomes  the  focus  of  a 
developing  TC.  Experiments  1 1  and  13  actually  developed  a  small-scale  vortex  but  it  was 
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not  resilient  and  instead  merged  with  the  larger  scale  circulation.  These  cases  are 
therefore  designated  as  developing  along  pathway  one.  Two  experiments,  14  and  16,  first 
developed  into  TC’s  along  pathway  one,  but  subsequently  weakened  and  later  underwent 
a  second  genesis  along  pathway  two.  These  are  therefore  designated  1/2  in  the  table.  Also 
shown  for  the  pathway  two  cases  is  the  maximum  averaged  tangential  wind  at  the  lowest 
model  level  that  occurs  in  the  period  before  the  SSCV  forms.  The  next  columns  give  the 
355  time  at  which  the  TC  develops  into  a  tropical  storm  (Tts)  using  the  criteria  that  the 
average  tangential  winds  at  the  lowest  model  level  reach  17.4  m  s'1,  the  RMW  at  this 
time,  the  time  at  which  the  TC  reaches  hurricane  strength  (Th)  using  a  criteria  of  33  m  s'1 
and  the  RMW  at  this  time. 

From  this  table  the  following  conclusions  can  be  drawn:  (1)  Both  pathways  are  about 
360  equally  represented  in  this  set  of  experiments.  (2)  The  moist  cases  develop  much  more 
rapidly  than  the  dry  cases  given  that  the  other  parameters  are  the  same.  (3)  Radiation 
tends  to  produce  faster  development.  (4)  Pathway  one  is  favored  if  the  initial  vortex  is 
strong  or  if  the  sea  surface  temperatures  are  lower.  (5)  A  significant  fraction  of  the 
systems  that  develop  along  pathway  two  reach  wind  speeds  of  approximately  10  m  s'1 
365  near  the  surface  prior  to  an  SSCV  forming,  so  they  are  already  a  borderline  tropical 
depression.  (5)  The  SSCV’s  that  form  for  systems  developing  along  pathway  two  have  a 
very  small  RMW  and  they  typically  develop  very  rapidly  into  tropical  storms  and 
hurricanes  undergoing  a  modest  growth  during  this  period.  On  the  other  hand,  systems 
developing  along  pathway  one  tend  to  have  a  larger  RMW  and  undergo  a  contraction  as 
370  they  develop.  Some  of  these  cases  develop  rapidly  also.  (6)  Not  unexpectedly,  larger 
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initial  vortices  tend  to  produce  larger  TCs.  (7)  Unexpectedly,  lower  SSTs  do  not  always 
result  in  a  slower  rate  of  development. 

For  systems  that  develop  along  pathway  two.  Table  3  shows  properties  of  the  “second 
mid-level  vortex”  at  the  time  of  formation  of  the  SSCV  and  the  maximum  wind  speeds 
375  and  radius  of  maximum  winds  at  the  surface.  For  these  RAMS  simulations,  we  refer  to 
the  strong  mid-level  vortex  that  forms  as  the  “second  mid-level  vortex”  since  the  original 
MCV  became  less  evident  as  the  low  level  winds  increased  during  the  early  evolution  of 
the  system.  The  second  mid-level  vortex  develops  winds  between  12-16  m  s'1  for  most 
cases  with  the  exception  of  Experiment  12  that  has  winds  of  19.1  m  s’1.  The  RMW  of  the 
380  mid-level  vortex  tends  to  be  relatively  large  for  the  dry  cases  and  the  large  initial  vortex 
cases,  ranging  from  47  to  91  km.  The  experiments  with  small  initial  vortices  and  lower 
SST’s  develop  notably  smaller  second  mid-level  vortices.  They  are  centered  at  a  height  of 
approximately  4. 3-5. 4  km.  Comparing  Table  2  with  Table  3  it  can  be  seen  that  the 
SSCV’s  attain  a  tangential  wind  speed  of  12  m  s’1  within  a  few  hours  of  forming.  Surface 
385  winds  tend  to  decrease  from  their  maximum  values  shown  in  Table  2  as  the  mid-level 
vortex  develops.  Surface  winds  at  the  time  the  SSCV  forms  for  the  dry  cases  are 
particularly  weak.  The  RMW  of  the  surface  winds  at  this  time  tends  to  be  significantly 
larger  than  the  RMW  of  the  mid-level  vortex.  It  is  particularly  large  for  the  dry  cases. 

390  3.2  Discussion  of  individual  cases 

Experiment  1 :  This  experiment  has  the  simplest  microphysics  with  only  the  cloud  water 
and  rain  categories  activated.  There  is  no  radiation  scheme  and  the  initial  vortex  is  small, 
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weak  and  has  a  moisture  anomaly.  Figure  2  a,  b  and  c  show  time  series  of  the  minimum 
surface  pressure,  the  maximum  azimuthally-averaged  tangential  wind  and  the  height  it 
occurs,  and  the  maximum  azimuthally  averaged  tangential  wind  at  the  model  level 
adjacent  to  the  surface  (z=98  m)  and  the  radius  it  occurs,  respectively.  The  minimum 
surface  pressure  shows  a  gradual  decrease  during  the  first  35  h  followed  by  a  substantial 
fall  during  the  next  15  h.  Genesis  occurs  at  approximately  t=39  h.  The  height  of  the 
400  maximum  tangential  winds  gradually  falls  after  12  h  becoming  close  to  the  surface  at  the 
time  of  genesis.  Near-surface  winds  started  to  increase  significantly  after  30  h.  There  is  a 
significant  contraction  of  the  radius  of  maximum  winds  after  30  h  to  approximately  15 
km  at  genesis.  It  remains  fairly  constant  during  the  next  12  h  while  the  winds  rapidly 
increase  to  hurricane  strength.  For  this  case  a  second  mid-level  vortex  did  not  form.  The 
405  system  clearly  developed  into  a  tropical  cyclone  along  pathway  one. 

Experiment  2:  This  case  is  identical  to  experiment  1  except  that  ice  physics  is 
included.  The  rate  of  development  was  considerably  slower.  Figure  3  shows  that  during 
the  first  twenty  four  hours  the  maximum  winds  shifted  to  near  the  surface  with  speeds 
reaching  in  excess  of  10  m  s’1.  A  second  prominent  midlevel  vortex  developed  after  60  h. 
410  The  sudden  formation  of  a  SSCV  with  a  radius  of  5  km  is  evident  in  Fig.  3c  at  t=80  h, 
that  was  concurrent  with  a  surface  pressure  fall  of  ~  2  mb.  Near-surface  tangential  winds 
of  the  SSCV  reach  12  m  s’1  a  couple  of  hours  later.  The  small  vortex  increased  in  size 
fairly  soon  after  forming  although  the  system  remained  relatively  small.  Figure  4a  and  b, 
show  vertical  cross  sections  of  the  averaged  tangential  component  of  the  horizontal 
415  velocity  at  t=79  and  83  h,  respectively.  Just  prior  to  the  formation  of  the  small  vortex  the 
midlevel  vortex  is  centered  at  approximately  5  km  above  the  surface  and  the  maximum 
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wind  speed  is  ~  12  ms’1  at  a  radius  of  ~35  km.  The  strongest  surface  winds  are  at  a  radius 
of  ~60  km.  A  few  hours  later  a  SSCV  has  formed  with  near  surface  wind  speeds  in  excess 
of  12  m  s'1.  During  this  period  and  afterwards  the  midlevel  vortex  continued  to  strengthen 
420  as  shown  in  Fig.  3b.  At  ~92  h  the  near  surface  winds  become  stronger  than  the  winds 
aloft.  By  this  time  they  are  already  at  tropical  storm  strength.  This  experiment  clearly 
demonstrates  the  importance  of  the  ice  phase  in  the  development  of  a  second  stronger 
mid-level  vortex  and  the  subsequent  formation  of  a  SSCV. 

Experiment  3:  This  case  is  the  same  as  Experiment  2  except  there  is  no  initial 
425  moisture  anomaly.  It  is  similar  to  an  experiment  conducted  by  N07  that  also  had  no 
moisture  anomaly  or  radiation.  In  this  experiment  the  environment  has  less  low  level 
moisture  than  used  in  the  experiments  conducted  by  N07  since  the  Jordan  sounding  was 
modified  as  discussed  in  section  2.2.  Figure  5  shows  that  the  rate  of  development  was 
extremely  slow  for  this  case.  A  SSCV  formed  at  t=147  h  and  a  pressure  fall  of  ~3  mb 
430  occurred  in  the  next  few  hours.  However,  the  fall  in  minimum  surface  pressure  was  only 
temporary.  Also,  the  near  surface  wind  speeds  decreased  to  less  than  10  m  s'1.  Sustained 
development  did  not  occur  until  after  160  h.  The  midlevel  vortex  was  considerably  larger 
than  for  Experiment  2  and  appears  to  be  related  to  the  more  widespread  distribution  of 
deep  convection  for  this  case,  whereas  the  presence  of  an  initial  vortex  moisture  anomaly 
435  for  Experiment  2  resulted  in  convection  being  more  focused  in  the  center.  Figure  6  shows 
the  average  relative  humidity  for  a  circular  area  with  radius  of  100  km  from  the  vortex 
center  at  z=2.6,  4.9  and  7.3  km.  The  vortex  moistens  at  low  levels  first  and  then  more 
gradually  at  upper  levels  due  to  vertical  transport  of  moisture  in  deep  convective  towers. 
At  approximately  t=110  h  all  three  levels  show  moistening  to  about  85%  of  saturation. 
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Subsequently,  the  moistening  increases  at  z=7.3  km,  remains  approximately  constant  at 
z=4.9  km  and  decreases  slightly  at  2.6  km.  It  appears  that  deep  moistening  was  an 
important  prerequisite  for  genesis,  however  it  is  interesting  that  genesis  did  not  occur 
until  a  considerable  period  of  time  after  high  values  of  relative  humidity  had  been 
reached  at  these  levels.  This  case  also  showed  less  spin  up  of  low-level  winds  prior  to  the 
445  development  of  a  second  midlevel  vortex  than  Experiment  2.  This  is  consistent  with  the 
results  of  N07  who  found  low-level  winds  did  not  increase  significantly  for  his 
experiments  that  did  not  include  an  initial  low-level  moisture  anomaly. 

Experiment  4:  This  case  is  similar  to  Experiment  2  except  that  the  size  of  the  initial 
vortex  and  moisture  anomaly  is  larger.  Development  was  similar  in  many  respects  to 
450  Experiment  2,  but  slightly  slower.  Additionally,  it  eventually  developed  a  significantly 
larger  tropical  storm  and  hurricane  than  for  Experiment  2.  Although  a  SSCV  formed,  it 
grew  substantially  in  size  as  the  system  intensified  (Table  2).  Why  this  system  developed 
into  a  larger  hurricane  than  other  systems  developing  along  pathway  two  is  not  certain, 
but  it  did  develop  quite  a  significant  large-scale  surface  circulation  prior  to  the  SSCV 
455  forming.  Convergence  of  this  large-scale  low-level  cyclonic  vorticity  may  have  been  a 
factor  influencing  the  ultimate  size  of  the  circulation  at  hurricane  strength. 

Experiment  5:  This  case  is  similar  to  Experiment  2  except  that  the  SST  is  reduced 
from  29  to  28  °C.  Figure  7  shows  genesis  occurring  at  approximately  t=95  h.  The 
pathway  to  genesis  is  different  than  for  Experiment  2.  A  significant  second  midlevel 
460  vortex  still  forms  although  it  is  smaller  since  early  convection  is  more  focused  near  the 
center  for  this  case  (not  shown).  There  is  a  more  gradual  decrease  in  the  radius  of 
maximum  winds  near  the  surface  rather  than  a  clear  jump  as  occurred  in  Experiment  2 
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when  the  very  small  vortex  formed.  When  the  radius  of  maximum  winds  near  the  surface 
reach  ~15  km,  wind  speeds  reach  tropical  cyclone  strength.  This  is  in  contrast  to 
465  Experiments  2,3  and  4,  which  show  the  sudden  formation  of  a  very  small  vortex  with  a 
radius  of  5-6  km  prior  to  genesis.  Although  genesis  occurs  along  pathway  one  for  this 
case  it  is  clearly  very  close  to  a  transition  point  whereby  a  small  increase  in  SST  would 
likely  result  in  development  along  pathway  two. 

Experiment  6:  This  case  is  a  larger  vortex  case  similar  to  Experiment  4  except  that  the 
470  SST  is  reduced  to  28  °C.  This  case  underwent  genesis  at  approximately  t=l  12  h  as  shown 
in  Figure  8.  At  this  time  a  very  small  region  of  strong  vorticity  developed  near  the  center 
accompanied  by  a  local  pressure  fall.  This  feature  evident  in  Figure  8c  was  not  persistent 
and  only  lasted  about  an  hour.  This  case  is  therefore  classified  as  developing  along 
pathway  one.  The  radius  of  maximum  winds  near  the  surface,  when  winds  reach  tropical 
475  cyclone  strength,  varies  considerably  with  an  average  of  ~40  km.  This  case  also 
developed  a  second  midlevel  vortex  prior  to  genesis.  This  case  is  similar  to  Experiment  5 
in  that  it  is  close  to  developing  along  pathway  two. 

Experiment  7:  This  case  is  similar  to  Experiment  5  except  that  the  SST  is  reduced 
further  to  27  °C.  Figure  9  shows  that  genesis  occurs  at  approximately  t=86  h.  This  case 
480  did  not  develop  a  second  midlevel  vortex  prior  to  genesis.  Later  as  the  system  developed 
into  a  tropical  storm  it  can  be  seen  in  Figure  9  b  that  mid-level  winds  were  often  stronger 
than  those  near  the  surface.  Even  though  the  system  developed  along  pathway  one  it  still 
formed  a  very  small  hurricane  comparable  to  many  of  those  that  developed  along 
pathway  two. 

485 
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3.3  Simulations  with  a  radiation  scheme 


In  this  section,  we  discuss  simulations  with  the  radiation  scheme  activated,  in  particular 
drawing  attention  to  differences  with  the  non-radiation  cases  discussed  previously. 
490  Experiment  8  is  similar  to  Experiment  2  except  that  the  Harrington  radiation  scheme  is 
included.  As  can  be  seen  in  table  2  this  led  to  a  significantly  faster  rate  of  genesis.  This  is 
true  in  general  for  simulations  that  include  radiation  compared  to  similar  simulations  that 
do  not.  The  effects  of  radiation  on  tropical  oceanic  systems  can  broadly  be  classified  as 
due  to:  (1)  Clear  sky  infrared  (IR)  cooling  (2)  Direct  radiation-convection  interaction,  and 
495  (3)  cloud-cloud  free  radiation  difference  (Gray  and  Jacobson,  1977).  Clear  sky  IR  cooling 

is  very  evident  in  the  environment  of  the  simulations  with  the  radiation  scheme  included. 
The  strongest  IR  cooling  occurs  at  9-10  km  above  the  surface,  which  is  consistent  with 
the  results  of  Zhang  and  Chou  (1999).  This  will  increase  the  convective  available 
potential  energy  (CAPE)  over  time  promoting  moist  convection.  This  is  probably  a 
500  significant  factor  responsible  for  the  faster  rate  of  genesis  for  the  cases  with  radiation.  IR 
cooling  causes  also  the  environmental  profile  to  become  closer  to  saturation,  which  may 
promote  also  the  earlier  development  of  convection.  Another  aspect  of  the  IR  cooling  is 
that  it  leads  to  an  increase  in  the  surface  pressure  in  the  environment  and  also  in  the  core 
of  the  system  early  on  when  latent  heating  due  to  deep  convection  is  relatively  small. 
505  This  is  particularly  evident  in  the  dry  simulations  since  it  takes  longer  for  deep 
convection  to  develop. 

Experiment  15  is  an  interesting  case  since  it  develops  along  pathway  two,  whereas 
Experiment  6,  which  is  identical  except  for  the  absence  of  radiation,  develops  along 
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pathway  one.  Figure  10  illustrates  this  development  and  includes  time  series  of  total 
liquid  and  ice  water  within  a  cylindrical  volume  with  a  radius  of  100  km  from  the  center, 
as  well  as  clear  sky  short  wave  radiation  at  the  surface.  For  the  radiation  case  the 
minimum  surface  pressure  increases  initially  due  to  IR  cooling  as  discussed  earlier,  in 
contrast  to  the  non-radiation  case  that  exhibits  a  decrease.  By  t=80  h  the  radiation  case 
has  developed  into  a  moderately  strong  hurricane,  whereas  the  case  without  radiation 
515  does  not  undergo  genesis  until  approximately  t=l  12  h  (Fig.  8).  Fig.  10  d  shows  that  the 
liquid  and  ice  content  is  far  larger  for  the  radiation  case.  There  are  large  pulses  in 
convective  activity  that  are  clearly  related  to  the  nighttime  absence  of  short  wave 
radiation.  A  modeling  study  of  the  effects  of  radiation  on  tropical  cyclones  by  Hobgood 
(1986)  suggests  that  the  diurnal  cycle  of  net  radiation  at  the  cloud  tops  steepens  the  lapse 
520  rate  at  night  and  increases  convection.  During  daylight  hours,  the  absorption  of  solar 
radiation  reduces  the  lapse  rate,  thus  reducing  convection.  This  direct  radiation- 
convection  interaction  may  be  responsible  for  the  diurnal  cycles  of  convective  activity 
that  occur  in  Experiment  15,  and  also  in  the  other  radiation  experiments  conducted  in  this 
study.  A  detailed  analysis  of  the  interaction  of  radiation  and  convection  is  beyond  the 
525  scope  of  this  study,  but  these  results  indicate  that  they  are  likely  to  have  an  important 
impact  on  the  rate  of  tropical  cyclogenesis  and  possibly  on  the  pathway  taken  for  this 
particular  posing  of  the  genesis  problem  within  an  MCV  embryo. 

The  second  mid-level  vortex  in  Experiment  15  develops  during  the  first  surge  in 
convective  activity  at  t=40  h.  This  second  mid-level  vortex  is  short  lived  relative  to  the 
530  non-radiation  case  (Fig.  10  d  and  Fig.  8  b).  A  SSCV  forms  only  eight  hours  after  the 
midlevel  vortex  develops,  at  approximately  9  am  local  time.  It  then  undergoes  sustained 
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development  through  the  daytime  hours.  Although  the  total  liquid  and  ice  in  the  system 
core  is  at  a  minimum  during  this  period  it  is  still  considerably  larger  than  the  non¬ 
radiation  case.  As  the  system  develops  the  radius  of  maximum  winds  of  the  small  vortex 
535  increases  from  ~5  km  at  t=49  h,  to  ~17  km  at  t=90  h  (Fig.  10c). 

The  three  radiation  simulations,  that  were  initialized  with  a  strong  initial  vortex, 
Experiments  10,  11,  and  13,  all  developed  quickly  along  pathway  one  and  formed 
relatively  large  tropical  cyclones.  Experiment  10  underwent  genesis  at  approximately 
t=30  h.  After  genesis  between  34-45  h,  the  midlevel  winds  were  generally  stronger  than 
540  at  low  levels  although  a  well-defined  secondary  vortex  did  not  form.  During  this  same 
period  the  near  surface  winds  increased  significantly  from  12  to  20  m  s’1.  Experiment  1 1, 
that  was  identical  to  Experiment  10  except  for  the  absence  of  an  initial  moisture  anomaly, 
underwent  genesis  at  approximately  t=45  h.  This  case  did  form  a  well-defined  second 
midlevel  vortex  after  genesis  at  t=47  h.  Also,  a  SSCV  formed  a  few  hours  later  at  t=52  h. 
545  This  was  concurrent  with  a  sudden  jump  in  the  radius  of  maximum  winds  that  decreased 
to  7  km.  However,  this  small  vortex  was  not  persistent  and  the  radius  of  maximum  winds 
increased  to  17  km  three  hours  later.  Therefore,  this  case  was  classified  as  developing  via 
pathway  one.  Experiment  13  underwent  genesis  at  approximately  t=35  h  evolving  along 
pathway  one.  Interestingly,  a  second  mid-level  vortex  developed  at  t=39  h  and  a  SSCV  at 
550  t=42  h.  This  small  vortex  was  weak  with  only  a  2  mb  decrease  in  surface  pressure  and 

there  was  not  a  concurrent  jump  in  the  radius  of  maximum  winds  when  it  formed.  It 
decayed  within  a  couple  of  hours. 

Experiments  14  and  16  are  similar  to  Experiments  5  and  7,  respectively,  except  that 
radiation  is  included.  These  lower  SST  cases  with  radiation  also  developed  along 
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555  pathway  one  initially,  but  subsequently  weakened  and  then  underwent  genesis  a  second 
time  along  pathway  two.  The  results  of  these  experiments  and  Experiment  15  discussed 
previously  suggest  that  inclusion  of  radiation  may  make  development  along  pathway  two 
more  likely.  However,  Experiment  17  that  uses  the  Chen-Cotton  radiation  scheme,  which 
is  identical  to  Experiment  8  otherwise,  showed  development  along  pathway  one.  The 
560  Chen-Cotton  scheme  is  simpler  than  Harrington’s,  having  less  realistic  microphysical- 
radiation  interactions.  Moreover,  comparing  the  clear  sky  IR  cooling  rates  showed  that 
they  were  significantly  larger  for  the  Chen-Cotton  scheme  than  for  the  Harrington 
scheme,  especially  at  lower  levels.  This  apparently  led  to  larger  values  of  CAPE 
developing  early  on,  as  well  as  bringing  the  air  closer  to  saturation  which  should  favor 
565  stronger  convective  activity  during  the  nighttime  between  t=9  to  21  hours.  Genesis 
occurred  during  this  first  pulse  in  convective  activity,  so  that  by  the  early  morning  the 
system  was  developing  rapidly  and  already  undergoing  intensification.  A  mid-level 
vortex  did  not  develop  for  this  case.  Although  results  for  the  Chen-Cotton  scheme  are  not 
likely  to  be  as  realistic  as  for  the  Harrington  scheme,  this  experiment  does  provide 
570  additional  evidence  that  radiation  may  have  significant  influence  on  the  pathway  taken  to 
genesis. 

3.4  The  relation  between  sea  surface  temperature  and  ice  production  on  the 
pathway  taken  to  genesis 

575 

Comparison  of  Experiments  1  and  2  clearly  demonstrate  the  importance  of  the  ice  phase 
in  promoting  development  along  pathway  two.  Furthermore,  experiments  show  that 
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lowering  SST  when  the  ice  phase  is  activated  can  cause  genesis  to  occur  along  pathway 
one  instead  of  two.  Development  along  pathway  two  is  favored  by  the  formation  of  a 
580  strong  mid-level  vortex,  which  is  related  to  a  substantial  stratiform  ice  layer  aloft  and  a 
prominent  mid-level  inflow.  Figure  1 1  shows  the  mixing  ratio  of  ice  and  radial  winds  for 
Experiment  two  at  75  h.  Between  30  to  40  km  from  the  center  there  are  high  ice 
concentrations  that  were  produced  in  deep  convective  towers.  The  base  of  the  ice  layer 
slopes  upwards  with  distance  form  the  center.  Beneath  this  base,  a  mid-level  inflow 
585  occurs  with  wind  speeds  of  approximately  2ms’1  at  a  height  of  7  km.  The  Coriolis  force 
acting  on  this  inflow  air  appears  to  be  mainly  responsible  for  the  spin-up  of  the  mid-level 
circulation  and  the  strong  vortex  evident  in  Figure  4. 

Given  the  importance  of  a  substantial  ice  layer  for  driving  a  mid-level  inflow  that 
produces  a  strong  mid-level  vortex  it  is  a  reasonable  hypothesis  that  lower  SST’s,  that  are 
590  likely  to  produce  weaker  convection  and  lesser  amounts  of  ice,  are  more  likely  to  favor 
development  along  pathway  one.  This  hypothesis  is  further  supported,  by  comparing  the 
amount  of  ice  and  vertical  mass  transport  for  Experiment  2,  which  has  an  SST  of  29  °C 
and  Experiment  5  that  is  identical  except  the  SST  is  28  °C.  Figures  12  a,  b  and  c  compare 
the  total  amount  of  ice  in  the  fine  grid  domain,  the  total  amount  of  ice  for  radius  less  than 
595  100  km  and  the  total  vertical  mass  transport  of  air  at  the  7.35  km  level  for  radius  less  than 

100  km.  For  Experiment  2,  that  has  a  higher  SST,  there  is  considerably  more  ice  in  the 
fine  grid  domain.  Also,  the  total  amount  of  ice  occurring  in  a  radius  less  than  100  km  is 
more  for  Experiment  2  although  it  is  not  so  notable.  There  are  large  oscillations  evident 
in  the  ice  content  particularly  for  Experiment  2  that  have  periods  of  approximately  ten 
600  hours.  The  total  vertical  mass  flux  at  the  7.35  km  level  for  radius  less  than  100  km  also 
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tends  to  be  higher  for  Experiment  2,  particularly  early  on.  Again  very  large  oscillations 
are  evident  that  are  discussed  further  in  section  3f.  These  results  indicate,  as  would  be 
expected,  that  higher  SST’s  tend  to  cause  stronger  convective  cells  that  produce  more  ice. 
The  outflow  layer  aloft  seen  in  Figure  lib  advects  ice  away  from  the  center  of  the 
605  system.  As  the  ice  falls  into  drier  environmental  air  at  mid  levels,  sublimation,  melting 
and  evaporation  result  in  cooling  that  would  contribute  to  driving  a  mid-level  inflow,  as 
will  be  discussed  in  section  4.  After  a  surge  in  convective  activity,  a  vertical  heating 
profile  more  similar  to  the  stratiform  region  of  MCSs,  which  is  characterized  by  latent 
heating  in  the  upper  troposphere  and  cooling  in  the  lower  troposphere  (Houze  1982; 
610  Johnson  and  Young  1983),  is  likely  to  become  more  prevalent,  even  in  the  center  of  the 
system.  A  linear  hydrostatic  analytic  model  of  the  forced  gravity  wave  response  to  a 
combination  of  both  convective  and  stratiform  heating  profiles,  was  derived  by  Nicholls 
et  al.  (1991).  The  stratiform  gravity  wave  mode  resulted  in  a  midlevel  inflow  and  lower- 
and  upper  level  outflows,  even  in  the  presence  of  the  convective  mode.  The  horizontal 
6 1 5  scale  of  these  circulations  expanded  away  from  the  heat  source  at  a  speed  depending  on 
the  Brunt- Vaisalla  frequency  and  the  vertical  wave  number  of  the  heating,  which  for  this 
case  would  be  approximately  20  m  s'1. 

3.5  Formation  of  the  small  surface  concentrated  vortex  for  cases  developing  along 
620  pathway  two 

Nine  of  the  experiments  listed  in  Table  2  underwent  genesis  along  pathway  two,  albeit 
two  of  these  cases  first  developed  along  pathway  one  before  weakening  and  then 


27 


625 


strengthening  again.  In  this  section  we  examine  the  conditions  existing  prior  to  the 
formation  of  the  SSCV  and  some  of  the  prominent  factors  that  appear  important  for  this 
process.  A  detailed  diagnosis  of  each  of  these  simulations  that  produced  the  SSCV  would 
be  desirable,  but  is  beyond  the  scope  of  this  exploratory  study. 

Experiment  2  formed  an  SSCV  at  t=80  h  with  winds  near  the  surface  reaching  12  m  s’1 
two  hours  later.  At  t=79  h  there  was  a  prominent  mid-level  vortex  (Fig,  4a)  and  high 
630  relative  humidity  throughout  the  troposphere  in  the  core  of  the  system.  These  conditions 
were  identified  by  N07,  as  important  precursors  for  the  formation  of  the  smaller-scale 
vortex.  Figure  13  shows  the  azimuthally-averaged  potential  temperature  perturbation  at 
t=79  h.  This  shows  a  significant  warm  anomaly  at  upper  levels  in  the  core  of  the  system. 
Beneath  this  warm  anomaly  there  is  a  cooling  between  z=2-4.5  km.  It  is  possible  that  this 
635  cooling  by  decreasing  the  low-level  stability  aided  the  development  of  convection  at  the 
center  a  short  time  later.  It  can  also  be  seen  that  there  is  a  mean  cooling  of  the  air  near  the 
surface  at  the  center,  which  would  not  be  conducive  to  the  development  of  a  convective 
cell  if  it  was  to  ingest  this  air.  During  the  next  hour  a  convective  cell  occurred  at  the 
center  and  an  SSCV  formed.  A  narrow  column  of  warm  air  extended  through  the  cold 
640  layer  and  the  surface  pressure  at  the  center  decreased  locally  by  ~  3  mb. 

The  situation  preceding  the  formation  of  the  SSCV  was  quite  complex.  Seven  hours 
earlier,  at  t=73  h,  there  was  a  peak  of  convective  activity  that  is  evident  in  the  time  series 
of  ice  content  and  vertical  mass  flux  shown  in  Fig.  12.  At  this  time  there  was 
considerable  convective  activity  southwest  of  the  center  that  created  strong  outflows  and 
645  vorticity  anomalies.  By  t=74  h  a  broad  region  of  cold  and  dry  air  lay  to  the  south  of  the 
center,  that  was  advecting  with  the  cyclonic  flow.  A  strong  cell  had  formed  at  the  eastern 
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edge  of  this  outflow  15  km  southeast  of  the  center  and  had  formed  a  local  region  of 
significantly  colder  and  drier  air  at  the  surface.  Convection  continued  on  the  western  edge 
of  the  broad  cold  and  dry  region  and  at  t=75  h  a  strong  cell  had  formed  15  km  due  south 
650  of  the  center.  There  were  also  active  cells  50  km  to  the  south  of  the  center  and  cells  30 
km  southwest  of  the  center.  By  this  time  the  very  cold  air  from  the  earlier  southeast  cell 
had  advected  nearer  the  center,  and  the  broad  region  of  cool  and  dry  air  at  the  surface  was 
in  the  southeast  quadrant.  At  t=76  h  cold  downdraft  air  from  the  cell  that  had  been  15  km 
south  of  the  center  was  contributing  to  form  a  pool  of  significantly  colder  and  drier  air 
655  extending  from  the  center  to  the  southeast  within  the  existing  broad  cold  air  region.  Also, 
at  this  time  new  development  had  occurred  ~5  km  south  east  of  the  location  of  the  earlier 
cell,  which  was  very  active  and  also  contributed  to  forming  the  cold  air  anomaly  in  this 
region.  There  was  significant  low-level  positive  vorticity  in  a  broad  region  southwest  of 
the  center,  presumably  a  result  of  the  numerous  convective  cells  that  had  occurred  earlier. 
660  There  was  a  local  surface  high  pressure  anomaly  beneath  the  new  active  cell,  but 
immediately  to  the  west  of  it  in  the  center  of  the  low  level  positive  vorticity  anomaly  a 
surface  mesolow  had  formed.  At  t=77h  the  mesolow  had  broadened  and  was  trailing  the 
surface  cold  region,  about  20  km  southeast  of  the  center.  Convection  in  this  region  and 
throughout  the  system  had  declined  by  this  time,  as  can  be  seen  in  Fig.  12.  This  decline  in 
665  convective  activity  was  probably  primarily  due  to  decreased  Convective  Available 
Potential  Energy  (CAPE),  as  the  low  level  air  was  cooled  and  dried  by  convective 
downdrafts.  The  region  of  significantly  colder  air  at  the  surface  that  was  ahead  of  the  low 
pressure  region  had  moved  closer  to  the  center  of  the  system  and  was  merging  with  the 
cold  air  that  had  previously  moved  to  the  center  from  the  earlier  southeast  cell  present  at 
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670  t=74  h.  So  the  center  had  become  significantly  cooler  than  a  few  hours  previously,  and  is 

the  main  reason  for  the  near  surface  cool  air  anomaly  seen  in  Fig.  13  at  t=79  h. 

Figure  14  a,  b,  c  and  d,  shows  horizontal  sections  of  vertical  vorticity,  potential 
temperature,  pressure,  at  z=98  m,  and  potential  temperature,  at  z=7.35  km,  for  t=  78  h, 
respectively.  The  vorticity  field  shows  a  broad  region  of  mainly  positive  vorticity 
675  anomalies  in  the  core  of  the  system.  Numerous  vorticity  filaments  can  be  identified.  They 
often  occur  at  boundaries  between  warm  and  cold  air.  The  relatively  cold  air  is  primarily 
caused  by  convective-scale  downdrafts.  The  low-level  air  near  the  center  of  the  system 
tends  to  be  colder  than  the  air  on  the  periphery.  The  mesolow  is  east  of  the  center  and  is 
located  on  the  edge  of  the  surface  cold  air  anomaly.  It  is  a  region  of  enhanced  positive 
680  vorticity.  Also  notable,  is  that  it  is  trailed  by  a  long  broad  filament  of  positive  vorticity 
that  circles  around  the  edge  of  the  central  cold  pool.  Fig.  14d  shows  that  the  warm  core 
aloft  is  asymmetrical  at  this  time.  The  position  of  the  warm  air  aloft  does  not  correlate 
well  with  the  mesolow  since  the  latter  is  a  shallow  feature.  There  is  a  warm  air  anomaly 
at  low  levels  within  the  mesolow  of  ~1  K  of  between  z=l  to  2  km.  At  this  time  the 
685  system  is  still  in  a  state  of  lowered  convective  activity  as  can  be  seen  in  Fig.  12.  There  are 
a  few  deep  cells  evident  by  the  high  values  of  potential  temperature  aloft  in  Fig.  14d,  at 
approximately  75  km  from  the  center,  but  very  little  activity  in  the  system  core.  During 
the  next  hour,  the  mesolow  broadened  and  moved  cyclonically  to  the  northeast, 
approximately  15  km  from  the  center,  remaining  on  the  edge  of  the  coldest  air  near  the 
690  center.  The  quiet  period  of  convective  activity  had  allowed  the  water  vapor  mixing  ratio 
near  the  surface  to  increase  to  relatively  high  values  outside  the  central  cold  pool. 
Consequently,  convective  activity  in  the  system  core  was  starting  to  become  significant 
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again.  At  this  time  there  was  a  surface  speed  maxima  approximately  10  km  due  north  of 
the  center  located  on  the  western  side  of  the  mesolow,  with  winds  from  the  north  having 
695  speeds  of  ~10  m  s'1  directed  towards  the  central  cold  pool.  It  was  at  this  location  on  the 
edge  of  the  central  cold  pool  that  the  cell  developed  that  formed  the  SSCV  an  hour  later. 

Figure  15  a,  b,  c,  d,  and  e,  shows  horizontal  sections  of  vertical  vorticity,  potential 
temperature,  water  vapor  mixing  ratio,  pressure,  at  z=98  m,  and  vertical  velocity,  at  z=1.6 
km,  respectively,  for  t=80  h.  The  SSCV,  which  has  just  formed  is  at  the  center  where 
700  there  is  a  very  small  negative  pressure  anomaly.  It  is  a  region  of  strong  positive  vorticity 
and  there  is  a  moderately  strong  updraft.  The  cell  ingested  relatively  warm  and  moist  air 
from  the  northwest.  The  cold  pool  of  air  adjacent  to  the  cell  is  responsible  for  the  near¬ 
surface  central  cold  anomaly  seen  in  the  azimuthally-averaged  potential  temperature  field 
shown  in  Fig.  13.  Low-level  vorticity  spiraled  into  the  cell  from  the  northwest  during  the 
705  next  hour,  which  appeared  to  aid  the  intensification  of  the  SSCV. 

For  Experiment  3,  the  origins  of  the  SSCV  could  be  traced  back  to  the  outflow  from  a 
strong  convective  cell  centered  70  km  to  the  east  of  the  center,  that  collided  with  a  pool 
of  cold  air  just  to  the  east  of  the  center.  This  caused  strong  easterlies  just  to  the  east  of  the 
center,  whereas  just  south  of  this  were  the  southwesterlies  of  the  cyclonic  circulation.  A 
710  VHT  with  a  weak  surface  low  formed  in  this  region  of  positive  low-level  vertical 
vorticity,  1 5  km  east  of  the  center.  It  remained  on  the  edge  of  the  more  central  cold  pool 
circling  once  around  the  center  and  within  three  hours  had  become  an  SSCV  that  was 
intensifying  rapidly.  By  t=151  h  tangential  wind  speeds  of  the  SSCV  had  reached  12  m  s’ 
,  but  a  cluster  of  cells  to  the  north  produced  a  strong  outflow  that  shortly  afterwards 
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715  swept  across  the  center  temporarily  weakening  the  SSCV,  as  can  be  seen  in  the  minimum 
surface  pressure  time  series  shown  in  Fig.  5. 

For  experiment  4,  a  cluster  of  cells  approximately  20  km  northwest  of  the  center 
produced  a  surface  cold  pool,  and  a  new  cell  formed  on  its  southern  edge  where  there  was 
significant  positive  vertical  vorticity.  This  VHT  had  a  small  surface  low  and  during  the 
720  next  two  hours  moved  cyclonically  to  15  km  southeast  of  the  center  intensifying  into  an 
SSCV.  During  this  period  a  squall  line  formed  west  of  the  center  that  rapidly  propagated 
southwards.  The  squall  line  did  not  impede  the  development  of  the  SSCV,  which  fed  on  a 
narrow  band  of  warm  and  moist  air  that  was  squeezed  between  the  older  cold  pool  to  the 
north  and  the  squall  line  outflow  to  the  south.  It  is  possible  that  the  strong  southwesterly 
725  surface  winds  just  south  of  the  cell  caused  by  the  squall  line  outflow  may  have 
contributed  to  the  vorticty  of  the  SSCV. 

For  Experiment  8,  a  strong  isolated  cell  20  km  from  the  center  created  a  cold  pool  of 
air  at  the  surface.  The  cyclonic  winds  that  were  stronger  aloft  due  to  the  strong  mid-level 
vortex  tilted  the  cell  down  shear,  so  that  the  rain  shaft  and  low  level  downdraft  was 
730  downstream  of  the  updraft.  At  this  time  the  outflow  winds  were  strong  and  the  cold  and 
dry  air  spread  out  at  the  surface  downstream  of  the  cell.  A  small  and  weak  surface  low 
formed  just  upstream  of  the  rain  shaft  and  it  was  a  region  of  notably  stronger  positive 
vertical  vorticity.  The  cell  trailed  the  cold  pool  around  the  center  for  the  next  three  hours 
feeding  on  low-level  air  from  its  right  rear  quadrant.  This  air  de-accelerated  after  entering 
735  the  low-pressure  region  and  fed  the  updraft  that  was  just  downstream  of  the  low.  As  the 
cold  pool  moved  closer  to  the  center  the  cell  moved  cyclonically  around  it  and  intensified 
to  form  an  SSCV  10  km  from  the  center. 
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For  Experiment  9,  the  evolution  was  significantly  different  than  the  other  experiments 
with  a  relatively  large  SSCV  forming  having  a  radius  of  1 1  km.  A  cluster  of  cells  80  km 
740  to  the  southwest  of  the  center  produced  a  broad  region  of  cold  outflow  with  strong 
southerly  winds.  A  weak  surface  mesolow  had  developed  west  of  the  center  with 
southerly  flow  on  its  eastern  side  from  the  cold  pool  outflow  and  northerly  flow  on  its 
western  side  from  the  large  scale  cyclonic  circulation.  Subsequently,  a  large  cluster  with 
several  cells  developed  explosively  in  this  region  50  km  west  of  the  domain  center.  Three 
745  of  the  more  intense  cells  that  had  formed  within  this  region  of  pre-existing  low-level 
vorticity  were  VHT’s  with  significant  cyclonic  circulation.  An  hour  later,  a  raised 
mesoscale  low  pressure  anomaly  with  a  cyclonic  circulation  had  been  created  between 
z=2-6  km.  In  the  decaying  phase  of  the  short-lived  MCS  the  raised  vortex  became 
stronger,  which  may  have  been  associated  with  the  development  of  local  mid-level  inflow 
750  as  low  level  cooling  due  to  melting  and  evaporation  of  hydrometeors  became  more 
prominent.  At  the  surface  a  strong  cold  pool  and  a  high  pressure  anomaly  formed  that 
generated  strong  outflow.  The  raised  mesoscale  vortex  was  a  persistent  feature  that 
moved  cyclonically  and  spiraled  towards  the  center  in  the  next  few  hours.  There  was  still 
a  significant  circulation  present  between  z=2-4  km  when  it  was  located  15  km  southwest 
755  of  the  center.  At  this  time  a  long  line  of  deep  convective  cells  had  developed  stretching 
from  west  to  east  through  the  center  of  the  fine  grid  domain.  Subsequently,  the  western 
line  moved  cyclonically  to  the  southwest,  and  the  eastern  part  of  the  line  to  the  northeast 
while  the  mesoscale  vortex  propagated  to  the  center.  Then  a  cluster  of  convection  at  the 
center  produced  cooler  air  and  a  higher  pressure  at  the  surface  while  intensifying  the 
760  mesoscale  vortex  aloft.  An  hour  later,  a  large  convective  cluster  explosively  developed  30 
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km  to  the  west  that  appeared  to  be  associated  with  the  earlier  line  of  convection  originally 
to  the  east  that  had  moved  cyclonically  round  the  center.  Low-level  convergence  at  the 
edge  of  the  spreading  central  cold  air  seems  to  have  also  promoted  its  development.  Two 
hours  later,  a  surface  low  had  formed  at  the  center.  Outflow  from  the  convection  to  the 
765  west  was  producing  strong  westerlies  on  its  southern  edge  that  may  have  aided  the 
development  of  the  cyclonic  winds  of  the  SSCV.  The  outflow  also  appeared  responsible 
for  initiating  a  convective  line  shortly  afterwards  just  south  of  the  center  that  rapidly 
propagated  to  the  southeast.  The  outflow  behind  the  convective  line  created  strong 
southwesterlies  on  the  south  side  of  the  SSCV  that  again  may  have  aided  its 
770  development.  Transient  multi-cellular  convection  occurred  around  the  edge  of  the  surface 
mesolow,  rather  than  near  its  center  as  the  SSCV  intensified.  It  is  possible  that  because 
there  was  a  significant  circulation  present  between  z=2-4  km  when  the  SSCV  developed 
that  convective  downdrafts  may  have  brought  some  of  the  cyclonic  vorticity  down  to  the 
surface.  There  were  however  significant  pockets  of  low-level  positive  vertical  vorticity 
775  near  the  surface  prior  to  the  SSCV  forming,  as  well  as  the  convective  outflows  that 
appear  to  favor  its  development,  so  it  probably  formed  mainly  from  the  bottom  up. 

For  Experiment  12  the  mid-level  vortex  was  very  strong  prior  to  the  SSCV  forming 
with  maximum  winds  of  19.1  m  s'1  and  a  large  radius  of  71  km.  The  origins  of  the  SSCV 
could  be  traced  back  to  a  small  weak  surface  low  with  a  positive  vertical  vorticity 
780  anomaly  25  km  southwest  of  the  center.  The  low  was  adjacent  to  a  local  cold  pool  and  a 
large  convective  cell  to  its  east.  Rainfall  from  the  convective  cell  was  producing  the  cold 
pool.  During  the  next  two  hours  the  small  low  moved  cyclonically  round  the  center 
trailing  the  cold  pool  showing  similarities  to  Experiment  8.  At  t=53  h  it  was  15  km  east 
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of  the  center.  There  was  no  significant  convection  present  in  this  region,  however  a 
strong  convective  cell  had  developed  65  km  to  the  northeast  that  was  producing  a  strong 
outflow.  At  t=  54  h,  the  pressure  in  the  small  low  had  fallen  significantly  and  it  had 
become  an  SSCV  25  km  north  of  the  center.  A  strong  convective  cell  was  adjacent  to  the 
low  level  vorticity  center  on  its  north  west  side.  Outflow  from  the  strong  cell  that  had 
been  to  the  northeast  was  reaching  this  location  and  enhancing  the  easterly  flow  on  the 
790  north  side  of  the  SSCV.  There  was  also  significant  positive  vorticity  anomalies  on  the 
outflow  boundary  of  the  cold  pool  that  appeared  to  be  merging  into  the  SSCV,  and  may 
have  contributed  to  its  development. 

For  Experiment  14  genesis  first  occurred  along  pathway  one  at  t=37  h,  followed  by 
weakening,  and  then  redevelopment  along  pathway  two,  forming  an  SSCV  at  t=66  h.  At 
795  t=57  h  the  cyclonic  circulation  at  the  surface  was  still  quite  strong  with  speeds  of  ~7  m  s'1 

between  a  radius  of  25-70  km  form  the  center.  A  mid-level  vortex  was  present  with  winds 
of  12  m  s'1  at  a  radius  of  35  km  and  a  height  of  6  km.  The  central  air  at  the  surface  was 
not  very  cold  at  this  time  and  the  water  vapor  mixing  ratio  was  high.  This  promoted 
convective  activity  during  the  next  few  hours  that  produced  a  broad  region  of  cold 
800  surface  air  at  the  center  of  the  system.  Although  convection  was  widespread  a  significant 
amount  was  close  to  the  center  and  this  appeared  to  be  responsible  for  increased  low  level 
convergence  and  spin  up  of  the  low  level  winds.  Also,  a  considerable  amount  of  ice  was 
produced  aloft  and  a  notable  mid-level  inflow  developed.  The  net  result  was  that  at  t=63 
h,  there  was  a  considerably  stronger  surface  circulation  with  maximum  mean  wind  speeds 
805  of  ~~9  m  s'1  at  a  radius  of  15  km  and  a  smaller  more  intense  mid-level  vortex  aloft.  At 
t=65  h,  a  strong  convective  cell  developed  10  km  east  of  the  center.  A  downdraft 
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occurred  downstream  on  its  northern  side  producing  a  local  cold  pool  and  a  small  surface 
pressure  formed  on  its  southern  side.  Strong  cyclonic  winds  moved  the  small  scale 
vorticity  10  km  to  the  north  in  the  next  hour  trailing  the  cold  pool  that  was  to  its  west,  and 
810  it  intensified  into  an  SSCV.  It  was  located  just  inside  the  wind  speed  maxima  of  the 
larger  scale  circulation  that  had  a  radius  of  15  km.  At  t=67  h  the  small  surface  low 
pressure  anomaly  was  situated  5  km  to  the  south  west  of  the  center  within  the  larger  scale 
low  pressure  region  and  surface  winds  in  the  circulation  that  the  SSCV  was  embedded  in 
had  reached  genesis  strength.  During  the  next  few  hours,  convection  near  the  center 
815  converged  the  relatively  strong  low  level  circulation  to  form  a  small  and  intense  vorticity 
center.  Although  an  SSCV  could  be  identified  during  the  second  time  it  underwent 
genesis,  it  did  show  some  characteristics  of  development  along  pathway  one,  so  this 
could  be  considered  a  marginal  case. 

For  Experiment  15  a  second  mid-level  vortex  of  only  moderate  size  developed,  with  a 
820  radius  of  47  km,  even  though  the  initial  vortex  was  large.  This  appeared  to  be  because 
convection  tended  to  be  more  focused  in  the  core  of  the  system,  possibly  because  the 
lower  SST  of  28  °C  led  to  a  reduction  of  outer  convection.  At  t=46  h  a  very  active  cluster 
was  located  50  km  to  the  east  and  shortly  after  a  strong  cell  also  developed  10  km  north 
of  the  center.  A  convective  downdraft  from  this  cell  contributed  to  the  central  cold 
825  region.  At  t=48  h  the  outflow  from  the  eastern  cluster  that  was  cyclonically  moving  round 
the  center  produced  strong  easterlies  at  approximately  25  km  north  of  the  center  that 
collided  with  the  central  cold  pool.  At  this  intersection,  where  there  was  strong  surface 
horizontal  wind  shear  a  cell  formed  that  quickly  developed  a  small  and  intense  low  level 
positive  vertical  vorticity  anomaly,  with  a  significant  surface  pressure  fall.  By  t=49  h,  the 
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830  SSCV  had  moved  to  10  km  northwest  of  the  center.  It  was  located  on  the  western  edge  of 
the  central  cold  pool.  The  rain  falling  from  an  active  cell  that  was  leading  the  small 
vorticity  maxima  was  evaporating  at  low  levels  and  contributing  to  cooling  the  air 
downstream.  The  cell  was  feeding  on  low-level  moist  air  coming  from  the  northwest. 

For  Experiment  16,  a  small  mid-level  vortex  formed  with  a  radius  of  only  23  km. 
835  Surface  winds  were  fairly  strong,  near  the  center  of  the  system  when  the  mid-level  vortex 
formed  because  it  had  previously  undergone  genesis  along  pathway  one,  similarly  to 
Experiment  14.  At  t=66  h  there  was  a  significant  cold  pool  at  the  center  produced  by 
recent  convective  activity  in  this  region.  Surface  winds  were  relatively  light  within  the 
cold  pool,  but  adjacent  to  its  edge  at  a  radius  of  30  km  there  was  a  significant  cyclonic 
840  circulation  with  wind  speeds  of  ~9  m  s'1.  An  almost  circular  filament  of  positive  vertical 
vorticity  ran  around  the  edge  of  the  central  cold  pool  at  this  time.  There  was  active 
convection  within  the  vorticity  filament  that  intensified  during  the  next  hour,  causing  it  to 
become  more  asymmetrical.  At  t=69  h  a  strong  cell  had  formed  15  km  south  of  the  center 
and  a  small  surface  low  pressure  was  forming  just  to  the  west  of  it.  The  cell  moved 
845  cyclonically  trailed  by  the  small  low.  It  was  inside  the  larger  scale  low  level  wind  speed 
maxima  that  had  a  radius  of  ~30  km.  There  was  a  small  but  significant  cold  pool  at  the 
center.  The  small  low  propagated  quickly  around  the  center  and  at  t=71  h  was  10  km 
southwest  of  the  center,  with  the  cold  pool  to  its  northeast.  A  closed  small-scale 
circulation  had  formed  around  it  at  the  surface  that  included  the  air  in  the  cold  pool  that 
850  had  strong  cyclonic  winds  around  the  low  pressure  center.  At  this  stage  it  was  an  SSCCV 
with  considerably  stronger  winds  than  the  larger  scale  circulation  that  could  still  be 
identified  at  a  radius  of  ~30  km. 
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3.6  Convective-scale  downdrafts  and  oscillations  in  convective  activity 

855 

Outflows  from  convective-scale  downdrafts  had  a  major  influence  on  the  formation  of 
SSCV’s,  as  discussed  in  the  previous  section.  The  cooling  and  drying  by  these 
downdrafts  was  typically  about  -3  K  and  -3  g  kg'1,  respectively.  This  is  slightly  less  than 
observations  of  boundary  layer  modification  by  strong  tropical  squall  lines  off  the  coast 
860  of  west  Africa  (Zipser  1977;  Fitzgarrald  and  Garstang  1981;  Johnson  and  Nicholls  1983; 

Nicholls  and  Johnson  1984).  One  possible  reason  for  this  is  because  prior  to  genesis  there 
is  considerable  moistening  in  the  core  of  the  system  that  would  tend  to  reduce 
evaporative  cooling. 

Figure  16  shows  time  series  of  the  average  potential  temperature  and  water  vapor 
865  mixing  ratio  within  a  distance  of  100  km  of  the  center,  at  z=305  m,  for  Experiment  2. 
Comparing  with  Figure  12  that  shows  the  ice  content  for  this  case,  there  is  a  close 
correlation  especially  for  the  water  vapor  mixing  ratio.  For  instance,  there  are  minima  of 
the  water  vapor  mixing  ratio  at  t=47,  54,  66,  75  and  84  h,  and  maxima  of  the  ice  content 
at  t=48,  56,  66,  74,  and  86  h.  This  suggests  that  as  convective  downdrafts  dry  the 
870  boundary  layer  air,  convective  activity  diminishes,  leading  to  less  ice  production.  There  is 
a  tendency  for  the  potential  temperature  anomalies  to  lag  the  moisture  anomalies.  One 
factor  that  may  contribute  to  this  lag  is  that  the  core  of  the  convective  downdrafts  at  this 
level,  are  often  warmer  relative  to  the  ambient  air.  This  is  possibly  because  precipitation 
drag  is  playing  a  significant  role  in  forcing  the  downdraft  in  addition  to  evaporative 
875  cooling,  so  that  it  is  negatively  buoyant,  in  spite  of  being  slightly  warmer  than  the 
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surrounding  air.  This  was  much  less  notable  at  the  lowest  model  level,  at  z=98  m,  where 
convective  downdrafts  were  almost  always  cooler  than  the  surrounding  air.  Another 
factor  is  that  light  rainfall  can  cause  cooling  and  moistening  at  low  levels,  rather  than 
drying.  Yet  another  factor,  is  that  compensating  subsidence  associated  with  a 
880  convectively  forced  gravity  wave  can  cause  warming  and  drying  (Nicholls  et  al.  1991; 
Nicholls  and  Pielke  2000). 

All  of  the  cases  without  radiation  showed  significant  oscillations  in  convective 
activity,  although  some  were  much  more  regular  than  others.  Periods  varied,  but  were 
typically  between  10-15  h.  After  a  period  of  intense  convective  activity  the  boundary 
885  layer  would  gradually  recover  and  Convective  Available  Potential  Energy  (CAPE)  would 
increase.  Convection  would  often  begin  fairly  close  to  the  center.  Spreading  outflows 
from  this  convection  would  tend  to  trigger  new  convection  further  from  the  center,  and  so 
convective  activity  would  work  its  way  outwards,  sometimes  as  partial  rings,  or  bands  of 
convection,  until  CAPE  was  diminished  again  over  a  broad  area.  Most  of  the  cases  with 
890  radiation  activated  showed  a  very  clear  diurnal  cycle  of  convective  activity,  as  seen  in 
Figure  10  d  and  e,  for  Experiment  15.  Experiment  17,  on  the  other  hand,  only  showed  one 
distinct  diurnal  oscillation  at  approximately  t=70  h,  followed  by  a  steady  rise  in  total 
condensate.  A  few  of  the  cases  showed  small  oscillations  in  between  the  diurnal  cycle. 
Because  of  the  large  amplitude  diurnal  component  in  these  simulations,  the  boundary 
895  layer  recovery  is  generally  prolonged.  However,  occasionally  recovery  can  take  place 
quickly  enough  that  another  oscillation  can  occur  during  the  daytime. 

An  SSCV  usually  formed  when  convective  activity  in  the  system  core  was  increasing, 
or  near  a  peak.  Exceptions  are  Experiments  4  and  12  that  formed  an  SSCV  during  a 
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minimum  in  convective  activity.  In  all  the  cases  simulated,  once  the  SSCV  had  formed  it 
900  developed  tangential  wind  speeds  of  12  m  s'1  within  2-3  hours.  For  the  cases  with 
radiation  activated,  intensification  to  tropical  storm  strength  occurred  only  a  few  hours 
later,  regardless  of  the  time  of  day  that  the  SSCV  formed.  Therefore,  development  of  the 
SSCV  did  not  appear  to  be  significantly  hindered  during  the  daytime  by  shortwave 
warming  of  the  upper  cloud  layer. 

905 

4  Discussion 

For  development  along  pathway  two,  the  fundamental  questions  are:  What  principally 
causes  the  strong  mid-level  vortex?  Is  it  important  to  the  formation  of  the  SSCV  and 
910  why?  These  questions,  cannot  be  definitively  answered  at  this  stage.  Nevertheless,  the 
experiments  that  were  conducted  do  shed  some  light  on  the  processes  that  are  likely  to  be 
relevant,  in  addition  to  those  proposed  by  N07.  Experiments  indicate  that  the  ice  phase  is 
important  for  producing  a  strong  mid-level  vortex  and  that  a  mid-level  inflow  at  the  base 
of  the  stratiform  ice  layer  precedes  its  development.  The  middle  level  inflow  could  be 
915  explained  by  axisymmetric  balance  dynamics  associated  with  the  radial  derivative  of  the 
latent  heating  due  to  a  “top-heavy”  heating  profile:  latent  heating  at  mid-  to  upper-  levels 
and  cooling  at  low  levels.  The  response  can  be  considered  a  slow  transverse  overturning 
circulation  to  the  heating  distribution.  For  a  12-h  time  average  M06  found  reasonable 
agreement  for  the  transverse  circulation  simulated  by  RAMS  and  Eliassen’s  balanced 
920  vortex  model  (Elissen  1951).  Alternatively,  since  the  mid-level  inflow  is  likely  to  be  the 
principal  cause  of  the  strong  mid-level  vortex,  and  during  its  development  there  will  be  a 
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certain  degree  of  imbalance,  another  perspective  is  provided  by  the  simple  linear 
hydrostatic  model  of  thermally-forced  gravity  waves,  that  neglects  the  Coriolis  force, 
developed  by  Nicholls  et  al.  (1991).  A  solution  was  obtained  for  a  combined  convective 
925  and  stratiform  heating  profile,  represented  by  a  half-sine  wave  between  the  surface  and 
the  tropopause  and  an  inverse-sine  wave,  respectively.  The  transient  response  consisted 
of  convective  and  stratiform  forced  gravity  wave  modes,  with  circulations  that  laterally 
expanded  away  from  the  heat  source  region,  at  a  speed  depending  on  the  Brunt  Vaisalla 
frequency  and  the  vertical  wave  number  of  the  heating.  For  the  slower  moving  stratiform 
930  mode  the  speed  would  be  approximately  20  m  s’1  for  the  tropical  atmosphere.  This  mode 
is  characterized  by  mid-level  inflow  and  lower-  and  upper-level  outflows.  Interestingly, 
even  though  the  combined  heating  profile  had  only  a  weak  net  low-level  cooling,  the 
stratiform  mode  that  was  superimposed  with  the  convective  mode  still  resulted  in 
significant  mid-level  inflow. 

935  Why  the  ice  phase  is  important  for  the  development  of  a  mid-level  inflow  and  why 
cases  that  produced  more  ice  appear  more  likely  to  evolve  along  pathway  two  remains  to 
be  answered.  The  liquid  phase  of  the  RAMS  model  includes  small  cloud  droplets  and 
rain.  There  is  no  sedimentation  of  cloud  droplets,  whereas  rain  falls  relatively  quickly.  On 
the  other  hand,  there  is  sedimentation  of  all  the  ice  categories  with  pristine  ice,  snow  and 
940  aggregates  generally  having  slower  fall  speeds  than  graupel  and  hail.  With  the  ice  phase 
included,  the  anvils  of  convective  towers  are  composed  mainly  of  ice,  rather  than  cloud 
droplets,  and  can  fall  to  lower  levels  and  undergo  sublimation  and  melting.  This 
difference  may  be  a  significant  factor  making  it  more  likely  that  a  stratiform  type  of 
heating  profile  develops.  Another  possibly  important  consideration  is  that  the  inflow  air 
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945  is  relatively  dry  since  it  descends  from  high  levels  and  also  because  some  of  it  originates 
from  the  periphery  of  the  system  that  has  not  been  moistened  so  much  by  vertical 
convective  moisture  transport  as  air  at  the  center.  This  would  tend  to  increase  sublimation 
and  evaporation  rates  of  hydrometeors  falling  into  the  inflow,  possibly  resulting  in  a 
positive  feedback.  Moreover,  the  inflow  advects  ice  at  the  base  of  the  stratiform  layer 
950  towards  the  center  of  the  system,  which  would  tend  to  focus  cooling  near  the  center, 
which  may  aid  in  the  contraction  and  spin  up  of  the  mid-level  vortex.  Another  major 
difference  when  ice  is  included  is  the  additional  latent  heat  of  fusion.  This  could  result  in 
enhanced  warming  aloft  and  cooling  at  low  levels,  increasing  the  tendency  for  mid-level 
convergence. 

955  Additional  factors  that  may  influence  the  spin-up  of  the  mid-level  vortex,  but  are 
probably  of  secondary  importance  to  stratiform  heating,  include:  (1)  Hydrometeor  drag 
that  could  enhance  mid-level  convergence  when  hydrometeors  are  concentrated  at  lower 
levels.  (2)  Mid-level  inflow  on  the  convective  scale  during  the  later  stage  of  convective 
cell  development.  At  this  time,  latent  heating  is  typically  skewed  towards  upper  levels,  as 
960  melting  and  evaporation  occur  at  low  levels.  The  development  of  a  convective  scale 
downdraft  at  low  levels  with  upward  motion  aloft  will  induce  mid-level  convergence. 
The  SSCV  typically  formed  near  to,  or  at  the  center  of  the  strong  mid-level  vortex.  If  it 
formed  away  from  the  center  it  would  move  cyclonically  and  spiral  into  the  center  of  the 
larger  scale  circulation,  usually  within  a  few  hours.  The  details  of  the  SSCV  formation 
965  showed  considerable  differences.  A  single  intense  VHT  often  played  a  role  in  its 
formation,  but  sometimes  multi-cellular  convection  was  present.  For  all  the  cases 
exhibiting  a  SSCV,  cold  pools  produced  by  convective  downdrafts  played  an  important 
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role.  Vertical  vorticity  was  typically  concentrated  at  the  edge  of  cold  pools  and  new 
convection  that  developed  in  these  regions  could  converge  and  stretch  this  pre-existing 
970  low-level  vorticity.  Also,  collisions  of  cold  pools  could  produce  convergence  and  regions 
of  enhanced  positive  vertical  vorticity  where  they  intersected  that  could  favor  the 
formation  of  an  SSCV.  A  notable  factor,  usually  present,  was  a  significant  cold  pool  at 
the  center  of  the  circulation  produced  by  previous  convective  activity.  An  SSCV  would 
often  form  along  its  boundary  where  vertical  vorticity  was  concentrated.  In  a  significant 
975  number  of  cases  a  precursor  to  a  SSCV  formed  approximately  20  km  from  the  center  of 
the  circulation,  when  a  strong  convective  cell  developed  above  a  region  of  pre-existing 
low-level  vertical  vorticity.  Because  of  the  stronger  cyclonic  circulation  of  the  mid-level 
vortex  aloft,  the  updraft  became  tilted  and  rain  falling  on  the  downstream  side  of  the  cell 
produced  a  cold  pool.  A  small  surface  low  pressure  formed  on  the  upstream  side  of  the 
980  cell  at  the  edge  of  the  cold  pool,  with  a  prominent  positive  vertical  vorticity  anomaly. 
This  could  strengthen  into  an  SSCV,  or  as  in  the  case  of  Experiment  2,  play  a  role  in  its 
formation  later  on.  The  small-scale  low  pressure  anomaly  would  trail  the  cold  pool, 
remaining  on  its  edge  as  it  spiraled  in  towards  the  center.  Intensification  of  convection 
would  often  occur  on  the  edge  of  the  pressure  low  as  it  reached  the  center  and  this  may 
985  have  been  facilitated  by  the  cold  air  at  mid-to-low  levels  at  the  center  of  the  mid-level 

vortex,  which  would  tend  to  decrease  low-level  static  stability.  The  convection  fed  on  a 
tongue  of  relatively  moist  and  warm  air  spiraling  towards  the  center  and  because  of  the 
reduced  low-level  static  stability  at  the  center  would  probably  be  in  a  region  of  increased 
low-level  CAPE.  The  SSCV  would  then  in  most  cases  develop  rapidly  into  a  tropical 
990  cyclone.  Therefore,  it  appears  that  the  formation  of  a  prominent  mid-level  vortex,  in 
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addition  to  causing  a  large  increase  in  the  efficiency  of  the  conversion  of  latent  heat 
energy  to  the  kinetic  energy  of  the  cyclonic  wind  field  as  found  by  N07,  may  play  a  role 
in  favoring  the  development  of  an  SSCV  due  to  its  vertical  wind  shear  and  decreased 
low-level  static  stability  at  its  center. 

995  A  complementary  interpretation  as  to  why  the  simulated  evolution  does  not 

continue  to  evolve  along  pathway  one  is  that  as  the  strong,  second  mid-level  vortex 
develops,  the  low  level  cooling  in  the  system  core  would  produce  a  divergent  tendency  at 
low  levels  that  might  arrest  the  contraction  of  the  surface  wind  speed  maxima.  The 
foregoing  thermodynamic  and  dynamic  interpretations  need  not  be  mutually  exclusive 
1000  and  future  work  should  quantify  their  relative  importance. 

Inclusion  of  active  radiation  is  found  to  significantly  increase  the  rate  of  genesis.  This 
appears  to  be  largely  explained  by  clear  sky  cooling  due  to  long  wave  radiation  at  upper 
levels  causing  CAPE  to  increase,  as  well  as  bringing  the  air  closer  to  saturation.  The 
Harrington  radiation  scheme  slightly  favored  evolution  along  pathway  two,  probably 
1005  because  of  cooling  aloft  promoting  ice  production;  whereas,  the  Chen-Cotton  scheme 
favored  evolution  along  pathway  one,  apparently  due  to  stronger  longwave  cooling, 
particularly  at  lower  levels,  causing  the  early  convective  activity  to  be  intense  enough  to 
quickly  spin  up  the  low-level  winds  to  tropical  cyclone  strength.  Additionally,  a  strong 
diurnal  oscillation  of  convective  activity  was  present  that  is  probably  mainly  due  to 
1010  daytime  warming  of  the  upper  cloud  layer  causing  an  increase  in  stability  aloft,  as  found 
by  Hobgood  (1986).  However,  even  in  the  complete  absence  of  radiation,  large  amplitude 
convective  oscillations  usually  occurred  with  a  shorter  period,  typically  between  10-15  h, 
that  were  associated  with  boundary  layer  modification  by  convective  downdrafts.  When 
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radiation  was  included  a  smaller  surge  in  convective  activity  sometimes  occurred  in  the 
1015  daytime,  between  diurnal  peaks.  This  indicates  that  the  boundary  layer  recovery  was 
often  significantly  shorter  than  the  diurnal  cycle.  Both  diurnal  and  semi-diurnal  cycles 
have  been  found  in  the  areal  extent  of  the  upper  level  cirrus  deck  of  fully  developed 
tropical  cyclones  (Browner  et  al.  1977;  Muramatsu  1983;  Lajoie  and  Butterworth  1984; 
Steranka  et  al.  1984;  Kossin  2002).  The  results  presented  here  for  the  early  stage  of  a 
1020  developing  system  may  be  relevant  to  these  previous  studies,  although  the  existence  of 
substantial  convective  downdrafts  and  cold  pools  in  the  core  of  developing  tropical 
cyclones  has  not  yet  been  established  definitively. 

5  Conclusions 

1025 

Results  of  idealized  numerical  simulations  suggest  that  there  are  two  canonical 
pathways  that  can  lead  to  tropical  cyclone  genesis  from  an  initial  weak  mid-level  vortex. 
One  pathway  proceeds  by  spin  up  of  surface  winds  that  remain  stronger  than,  or  of 
comparable  magnitude  to  winds  aloft  until  they  become  large  enough  that  they  reach 
1030  tropical  depression  strength,  defined  herein  as  wind  speeds  exceeding  12  m  s'1.  Vorticity 
is  concentrated  in  numerous  short-lived  rotating  deep  convective  towers  that  collectively 
drive  a  system  scale  inflow  at  low  levels.  Vorticity  gradually  builds  in  the  central  region 
as  the  system  scale  inflow  produces  increasing  cyclonic  vorticity  and  as  small-scale 
cyclonic  vorticity  anomalies  generated  by  the  VHTs  are  converged  at  low  levels  and 
1035  aggregate.  This  is  usually  a  stochastic  process  with  no  prominent  events  leading  to 
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sudden  falls  of  the  surface  pressure.  The  characteristics  of  pathway  one  are  consistent 
with  the  stochastic  view  of  genesis  as  articulated  by  Ooyama  (1982). 

The  second  pathway  at  early  times  proceeds  in  the  same  manner  as  the  first,  but 
before  surface  winds  reach  tropical  cyclone  strength  the  mid-level  winds  increase 
1040  significantly  and  become  stronger  than  the  near-surface  tangential  winds.  The  low-level 
air  becomes  relatively  cold  (between  approximately  1.5-5  km)  in  the  core  as  the  second 
stronger  mid-level  vortex  develops.  This  is  then  followed  by  the  sudden  formation  of  a 
smaller-scale  vortex  at  the  center  of  the  larger  scale  circulation.  The  small  vortex  is  warm 
core  and  produces  a  narrow  column  of  warm  air  within  the  relatively  colder  low-level  air 
1045  associated  with  the  mid-level  vortex.  It  extends  upwards  into  the  larger-scale  warm  core 
aloft.  Its  rapid  formation  often  leads  to  a  sudden  local  decrease  of  the  minimum  surface 
pressure  of  a  few  millibars.  Once  formed  it  is  resilient  and  becomes  the  focus  of  a 
strengthening  and  widening  vorticity  center.  Sometimes  the  pressure  continues  to  fall 
immediately  following  this  event  and  the  system  intensifies  rapidly.  Evolution  along 
1050  pathway  two  has  similarities  to  the  conceptual  model  of  Bister  and  Emanuel  (1997;  Fig. 
13,  section  8)  that  places  emphasis  on  the  role  of  the  mid-level  vortex.  However,  the 
development  of  low-level  vorticity  in  these  simulations  does  not  appear  to  be  primarily 
due  to  a  top  down  process,  as  they  envision.  The  surface  wind  speeds  of  some  of  the 
simulated  systems  evolving  along  pathway  approach  or  reach  tropical  cyclone  strength 
1055  prior  to  a  prominent  mid-level  vortex  developing  so  that  the  subsequent  formation  of  a 
SSCV  in  these  cases  could  be  considered  intensification  rather  than  genesis. 

Systems  that  evolved  along  pathway  two  developed  a  second  mid-level  vortex  that 
was  considerably  stronger  than  typically  occurs  in  the  stratiform  region  of  an  MCS,  with 
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wind  speeds  in  excess  of  12  m  s'1.  For  the  experiments  conducted,  its  horizontal  scale  as 
1060  measured  by  the  radius  of  maximum  tangential  winds  ranged  from  23  to  91  km,  at  the 
time  the  SSCV  formed.  It  was  shown  that  the  ice  phase  played  an  important  role  in 
producing  a  strong  second  mid-level  vortex.  Mid-level  inflow  that  developed  at  the  base 
of  the  ice  layer  appeared  to  be  the  main  factor  responsible  for  its  formation.  Results 
suggest  that  higher  SST’s  lead  to  more  ice  production,  which  favors  the  development  of  a 
1065  strong  mid-level  vortex  and  evolution  along  pathway  two. 

There  were  significant  differences  in  the  way  that  the  SSCV  formed  for  cases  that 
developed  along  pathway  two.  In  all  cases  there  appeared  to  be  significant  stretching  of 
pre-existing  low-level  vertical  vorticity  anomalies  by  convective  scale  updrafts.  Cold 
pools  produced  by  convective  scale  downdrafts  played  important  roles  in  the  formation 
1070  process.  In  some  cases,  the  vertical  wind  shear  associated  with  the  strong  midlevel  vortex 
influenced  the  convective  activity  and  development  of  a  small  low  level  vorticity 
anomaly  that  would  later  either  evolve  into  the  SSCV,  or  play  a  part  in  its  formation. 

Implementation  of  active  radiation  in  the  model  via  the  Harrington  scheme  increased 
the  tendency  to  develop  along  pathway  two.  In  contrast,  the  Chen-Cotton  radiation 
1075  scheme  increased  the  tendency  to  develop  along  pathway  one.  Both  schemes  led  to  an 
increase  in  the  rate  of  genesis  compared  to  non-radiation  cases  that  were  otherwise 
identical.  Strong  diurnal  oscillations  in  convective  activity  occurred  with  radiation 
activated.  SSCV’s  that  formed  during  the  daytime  when  net  convective  activity  was 
reduced  did  not  appear  to  be  significantly  hindered  from  developing. 

1080  The  RAMS  results  indicate  that  this  initial  value  problem  may  be  near  a  bifurcation 
point,  so  that  changes  to  model  physics  could  have  considerable  impact  on  the  outcome 
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of  experiments  and  the  likelihood  of  development  along  a  particular  pathway.  There  are 
outstanding  questions  regarding  the  physical  processes  involved  in  the  formation  of  the 
strong  mid-level  vortex  and  the  SSCV,  that  occur  during  development  along  pathway 
1085  two.  In  particular,  the  role  of  the  mid-level  vortex  in  the  sudden  formation  of  an  SSCV 
needs  to  be  further  understood.  At  present,  there  is  no  conclusive  observational  evidence 
of  a  tropical  cyclone  undergoing  genesis  along  pathway  two,  and  so  for  the  time  being  the 
existence  of  this  pathway  remains  hypothetical,  being  based  on  the  results  of  numerical 
modeling  studies  with  WRF  (N07)  and  with  RAMS  in  this  study.  The  uncertainty  in  the 
1090  realism  of  the  second  pathway  simulated  in  these  idealized  configurations  and  in  N07 
offer  some  motivation  to  examine  recent  observational  data  collected  during  the  Tropical 
Cyclone  Structure-2008  (TCS-08)  and  the  Pre-Depression  Investigation  of  Cloud 
Systems  in  the  Tropics  (PREDICT)  experiments  to  examine  if  some  TCs  do  indeed 
develop  in  the  manner  of  pathway  two. 

1095 

Ac  know!  edgem  ents: 

We  are  grateful  to  Saurabh  Barve  for  providing  computational  assistance,  and  Thomas 
Cram  and  Dr.  Wesley  Terwey  for  aid  with  data  analysis  code.  This  work  was  supported 
in  part  by  the  National  Science  Foundation,  under  Grants  NSF  AGS  0965721  and  ATM- 
1 100  0733380,  and  by  the  U.S.  Office  of  Naval  Research  Grant  No.  N00014-03-1-0185. 


48 


1105  References 


Bister,  M.,  and  Emanuel,  K.  A.:  The  genesis  of  hurricane  Guillermo:  TEXMEX 
analyses  and  a  modeling  study.  Mon.  Wea.  Rev.,  125,  2662-2682,  1997. 

Braun,  S.  A.:  A  cloud-resolving  simulation  of  Hurricane  Bob  (1991):  Storm 
1110  structure  and  eyewall  buoyancy.  Mon.  Wea.  Rev.,  130,  1573-1592,  2002. 

Braun,  S.  A.,  M.  T.  Montgomery,  K.  J.  Mallen,  and  Reasor,  P.  D.:  Simulation  and 
interpretation  of  the  genesis  of  tropical  stonn  GERT(2005)  as  part  of  the  NASA 
tropical  cloud  systems  and  processes  experiment.  J.  Atmos.  Sci.,  67,  999-1025,  2010. 

Browner,  S.  P.,  W.  L.  Woodley,  and  Griffith,  C.  G.:  Diurnal  oscillation  of  the  area 
1115  of  cloudiness  associated  with  tropical  storms.  Mon.  Wea.  Rev.,  105,  856-864,  1977. 

Charnock,  H.:  Wind  stress  on  a  water  surface.  Quart.  J.  Roy.  Met.  Soc.,  81,  639,  1955. 

Chen  C.,  and  Cotton,  W.  R.:  The  physics  of  the  marine  strato-cumulus-capped 
mixed  layer.  J.  Atmos.  Sci.,  44,  2951-2977,  1987. 

Clark,  T.L.,  and  Farley,  R.  D.:  Severe  downslope  windstorm  calculations  in  two  and 
1120  three  spatial  dimensions  using  anelastic  grid  nesting:  A  possible  mechanism  for 
gustiness.  J.  Atmos.  Sci.,  41,  329-350,  1984. 

Davis,  C.  A.,  and  Bosart,  L.  F.:  Numerical  simulations  of  the  genesis  of  Hurricane 
Diana  (1984).  Part  1:  Control  Simulation.  Mon.  Wea.  Rev.,  129,  1859-1881,  2001. 

Dunkerton,  T.  J.,  M.  T.  Montgomery,  and  Z.  Wang,  2009:  Tropical  cyclogenesis  in  a 

1 125  tropical  wave  critical  layer:  Easterly  waves.  Atmos.  Chem.  Phys.,  9,  5587-5646. 

Fitzjarrald,  D.  R.,  and  Garstang,  M.:  Vertical  structure  of  the  tropical  boundary  layer. 
Mon. Wea.  Rev.,  109, 1512-1526,  1981. 


49 


Gray,  W.  M.  and  Jacobson,  R.  W.  J.:  Diurnal  variation  of  deep  cumulus  convection. 
Mon.  Wea.  Rev.,  105,  1171-1188,  1977. 

1 130  Hack,  J.  J.,  and  Schubert,  W.  H.:  Nonlinear  response  of  atmospheric  vortices  to 

heating  by  organized  cumulus  convection.  J.  Atmos.  Sci.,  43,  1559-1573,  1986. 

Harr,  P.  A.,  and  Elsberry,  R.  L.:  structure  of  a  mesoscale  convective  system 

embedded  in  Typhoon  Robya  during  TCM-93.  Mon.  Wea.  Rev.,  124,  634-652,  1996. 

Harr,  P.  A.,  M.  S.  Kalafsky,  and  Elsberry,  R.  L.:  Environmental  conditions  prior  to 

1 135  formation  of  a  midget  tropical  cyclone  during  TCM-93.  Mon.  Wea.  Rev.,  124,  1693- 
1710,  1996. 

Harrington,  J.  Y.,  G.  Feingold,  and  Cotton,  W.  R.:  Radiative  impacts  on  the  growth  of  a 
population  of  drops  within  summertime  arctic  stratus.  J.  Atmos.  Sci.,  57,  766-785, 
2000. 

1 140  Hendricks  E.  A.,  M.  T.  Montgomery,  and  Davis,  C.  A.:  On  the  role  of  “vortical”  hot 
towers  on  tropical  cyclone  formation.  J.  Atmos.  Sci.,  61,  1200-1232,  2004. 

Hill,  G.E.:  Factors  controlling  the  size  and  spacing  of  cumulus  clouds  as  revealed 
by  numerical  experiments.  J.  Atmos.  Sci.,  31,  646-673,  1974. 

Hobgood,  J.  S.:  A  possible  mechanism  for  the  diurnal  oscillation  of  tropical 

1145  cyclones.  J.  Atmos.  Sci.,  43,  2901-2922,  1986. 

Houze,  R.  A.,  W-C  Lee,  and  Bell,  M.  M.:  Convective  contribution  to  the  genesis  of 
Hurricane  Ophelia.  Mon.  Wea.  Rev.,  137,  2778-2800,  2009. 

Johnson,  R.  H.  and  Nicholls,  M.  E.:  A  composite  analysis  of  the  boundary  layer 
accompanying  a  tropical  squall  line.  Mon.  Wea.  Rev.,  1 11,  308-319,  1983. 

1150  Jordan,  C.  L.:  Mean  soundings  for  the  West  Indies  area.  J.  Meteor.,  15,  91-97,  1958. 


50 


Klemp,  J.  B.,  and  Wilhelmson,  R.  B.:  The  simulation  of  three-dimensional  convective 
storm  dynamics.  J.  Atmos.  Sci.,  35,  1070-1086,  1978. 

Kossin,  J.  P.:  Daily  hurricane  variability  inferred  from  GOES  infrared  imagery. 

Mon.  Wea.  Rev.,  130,  2260-2270,  2002. 

1155  Lajoie,  F.  A.,  and  Butterworth,  I.  J.:  Oscillation  of  high-level  cirrus  and  heavy 

precipitation  around  Australian  region  tropical  cyclones.  Mon.  Wea.  Rev.,  1 12,  535- 
544,  1984. 

Lilly,  D.  K.:  On  the  numerical  simulation  of  buoyant  convection.  Tellus,  14,  148- 
172,  1962. 

1 160  Liu,  Y.,  D.-L.  Zhang,  and  Yau,  M.  K.:  A  multiscale  numerical  study  of  Hurricane 

Andrew  (1992).  Part  1:  Explicit  simulation  and  verification.  Mon.  Wea.  Rev.,  125, 
3073-3093,  1997. 

Louis,  J.-F.:  A  parametric  model  of  vertical  eddy  fluxes  in  the  atmosphere.  Bound. -Layer 
Meteor.,  17,  187-202,  1979. 

1 165  Mapes,  B.  E.,  and  Houze,  R.  A.:  Diabatic  divergence  profiles  in  western  Pacific 
mesoscale  convective  systems.  J.  Atmos.  Sci.,  52,  1807-1828,  1995. 

Meyers,  M.  P.,  R.  L.Walko,  J.  Y.  Harrington,  and  Cotton,  W.  R.:  New  RAMS  cloud 
microphysics  paramererization.  Part  II:  The  two-moment  scheme.  Atmos.  Res.,  45,  3- 
39,  1997. 

1 170  Montgomery,  M.  T.,  M.  E.  Nicholls,  T.  A.  Cram,  and  Saunders,  A.  B.:  A  vortical 
hot  tower  route  to  tropical  cyclogenesis.  J.  Atmos.  Sci.,  63,  355-386,  2006. 

Montgomery,  M.  T.,  Z.  Wang,  and  Dunkerton,  T.  J.:  Coarse,  intermediate  and  high 
resolution  numerical  simulations  of  the  transition  of  a  tropical  wave  critical  layer  to  a 


51 


tropical  storm.  Atmos.  Chem.  andPhys.,  10,  10803-10827,  2010. 

1 175  Montgomery  M.  T.,  and  Smith,  R.  K.:  Tropical  cyclone  formation:  Theory  and 

idealized  modeling:  Report  for  the  Seventh  International  Workshop  on  Tropical 
Cyclones,  La  Reunion,  Nov.  2010.  World  Meteorological  Organization,  Geneva, 
Switzerland,  2010. 

Montgomery,  M.  T.,  C.  Davis,  T.  Dunkerton,  Z.  Wang,  C.  Velden,  R.  Torn,  S.  J. 

1180  Majumdar,  F.  Zhang,  R.  K.  Smith,  L.  Bosart,  M.  M.  Bell,  J.  S.  Haase,  A.  Heymsfield, 

J.  Jensen,  T.  Campos  and  Boothe,  M.  A.:  The  Pre-Depression  Investigation  of 
Cloud  Systems  in  the  Tropics  (PREDICT)  Experiment:  Scientific  Basis,  New  Analysis 
Tools,  and  Some  First  Results.  Bull.  Amer.  Meteor.  Soc.,  93,  153-172,  2012. 

Montgomery  M.  T.,  and  Smith,  R.  K.:  The  genesis  of  Typhoon  Nuri  as  observed  during 
1185  the  Tropical  Cyclone  Structure  2008  (TCS08)  field  experiment.  Part  2:  Observations  of 
the  convective  environment.  Atmos.  Chem.  Phys.,  11,  31115-31136,  2011. 

Muramatsu,  T.:  Diurnal  variations  of  satellite-measured  Tbb  areal  distribution  and 
eye  diameter  of  mature  typhoons.  J.  Meteor.  Soc.  Japan,  61,  77-89,  1983. 

Nicholls,  M.  E.,  and  Johnson,  R.  H.:  A  model  of  a  tropical  squall  line  boundary  layer 
1 190  wake.  J.  Atmos.  Sci.,  41,  2774-2792,  1984. 

Nicholls,  M.E.,  and  R.A.  Pielke,  and  Cotton,  W.  R.:  Thermally  forced  gravity  waves 
in  an  atmosphere  at  rest.  J.  Atmos.  Sci.,  48,  1869-1884,  1991. 

Nicholls,  M.E.,  and  Pielke,  R.  A.:  Thermally-induced  compression  waves  and 

gravity  waves  generated  by  convective  storms.  J.  Atmos.  Sci.,  57,  3251-3271,  2000. 
1195  Nolan,  David  S.:  What  is  the  trigger  for  tropical  cyclogenesis?  Aust.  Met.  Mag.  56, 
241-266,  2007. 


52 


Ooyama  K.,  1969:  Numerical  simulation  of  the  life  cycle  of  tropical  cyclones.  J.  Atmos. 


Sci.,  26,  3-40,  1969. 

Pandya,  R.  E.,  and  Durran,  D.  R.:  The  influence  of  convectively  generated  thermal 

1200  forcing  on  the  mesoscale  circulation  around  squall  lines.  J.  Atmos.  Sci.,  53,  2924- 
2951,  1996. 

Powell  M.  D.:  Boundary-layer  structure  and  dynamics  in  outer  hurricane  rainbands.  Part 
II:  Downdraft  modification  and  mixed  layer  recovery.  Mon,  Wea.  Rev.,  118,  918-938, 
1990. 

1205  Raymond,  D.  J.,  C.  Lopez-Carrillo,  and  Lopez  Cavazos,  I.:  Case  studies  of 

developing  east  Pacific  easterly  waves.  Quart.  J.  Roy.  Meteor.  Soc.,  124,  2005-2034, 
1998. 

Raymond,  D.  J.,  and  Lopez-Carrillo,  C.:  the  vorticity  budget  of  Typhoon  Nuri  (2008). 
Atmos.  Chem.  Phys.,  10,  16589-16635,  2010. 

1210  Raymond,  D.  J.,  S.  L.  Sessions,  and  Lopez  Carrillo,  C..  Thermodynamics  of  tropical 
cyclogenesis  in  the  northwest  Pacific.  J.  Geophys.  Res.,  D18101,  18,  PP.,  2011 
doi:  10. 1029/201 1JD015624,  2011. 

Reasor,  P.,  M.  T.  Montgomery,  and  Bosart,  L.:  Mesoscale  observations  of  the 
genesis  of  hurricane  Dolly  (1996).  J.  Atmos.  Sci.,  62,  3151-3171,  2005. 

1215  Ritchie,  E.  A.,  and  G.  J.  Holland,  G.  J.:  Scale  interactions  during  the  formation  of 
Typhoon  Irving.  Mon.  Wea.  Rev.,  125,  1377-1396,  1997. 

Rogers,  R.  L.,  S.  S.  Chen,  J,  E,  Tenerelli,  and  Willoughby,  H.  E.:  A  numerical  study 
of  the  impact  of  vertical  shear  on  the  distribution  of  rainfall  in  Hurricane  Bonnie 
(1998).  Mon.  Wea.  Rev.,  131,  1577-1599,  2003. 


53 


1220  Schubert,  W.  H.,  J.  J.  Hack,  P.  L.  Silva  Dias,  and  Fulton,  S.  R.:  Geostrophic 
adjustment  in  an  axisymmetric  vortex.  J.  Atmos.  Sci.,  37,  1464-1484,  1980. 

Simpson,  J.,  E.  Ritchie,  G.  J.  Holland,  J.  Halverson,  and  Stewart,  S.:  Mesoscale 
interactions  in  tropical  cyclogenesis.  Mon.  Wea.  Rev.,  125,  2643-2661,  1997. 

Sippel,  J.  A,  J.  W.  Nielsen-Gammon,  Allen,  S.  E.:  The  multiple-vortex  nature  of 
1225  tropical  cyclogenesis.  Mon.  Wea.  Rev.,  134  (7),  1796-1814,  2006. 

Smagorinsky,  J.  S.:  General  circulation  experiments  with  the  primitive  equations. 

1:  The  basic  experiment.  Mon.  Wea.  Rev.,  91,  99-164,  1963. 

Smith,  R.  K.,  and  Montgomery,  M.  T.:  Observations  of  the  convective  environment 
in  developing  and  non-developing  tropical  disturbances.  Q.  J.  R.  Meteorol.  Soc.  In 
1230  review,  2012. 

Steranka,  J.,  E.  B.  Rodgers,  and  Gentry,  R.  C.:  The  diurnal  variation  of  Atlantic 

ocean  tropical  cyclone  cloud  distribution  inferred  from  geostationary  satellite  infrared 
measurements.  Mon.  Wea.  Rev.,  112,  2338-2344,  1984. 

Tripoli,  G.J.,  and  Cotton,  W.  R.:  The  use  of  ice-liquid  water  potential  temperature  as 
1235  a  thermodynamic  variable  in  deep  atmospheric  models.  Mon.  Wea.  Rev.,  109,  1094- 
1102,  1981. 

Walko,  R.  L,  W.  R.  Cotton,  J.  L.  Harrington,  Meyers,  M.  P.:  New  RAMS  cloud 

microphysics  parameterization.  Part  I:  The  single-moment  scheme.  Atmos.  Res.,  38, 
29-62,  1995. 

1240  Wang,  Z.,  M.  T.  Montgomery,  and  Dunkerton,  T.  J.:  Genesis  of  Pre-hurricane  Felix 
(2007).  Part  1:  The  role  of  the  easterly  wave  critical  layer.  J.  Atmos.  Sci.,  67,  1711- 
1729,  2010. 


54 


Zehr,  R.  M.:  Tropical  cyclogenesis  in  the  western  North  Pacific.  NOAA  Tech.  Rep. 
NESDIS  61,  181  pp.  [Available  from  U.S.  Department  of  Commerce, 

1245  NOAA/NESDIS,  5200  Auth  Rd.,  Washington,  DC  20233.],  1992. 

Zhang,  C,  and  Chou,  M-D.:  Variability  of  water  vapor,  infrared  radiative  cooling, 
and  atmospheric  instability  for  deep  convection  in  the  equatorial  western  Pacific.  J. 
Atmos.  Sci.,  56,  711-723,  1999. 

Zhang,  D.  -L.,  L.  Tian,  and  M.  -J.  Yang,  M.-J.:  Genesis  of  Typhoon  Nari  (2001)  from  a 
1250  mesoscale  convective  system.  J.  of  Geophys.  Res.,  doi:  10. 1029/201 1JD016640,  2011. 
Zipser,  E.  J.:  Mesoscale  and  convective-scale  downdrafts  as  distinct  components  of 
squall-line  structure.  Mon.  Wea  Rev.,  105,  1568-1589,  1977. 

Zipser  E.  J.,  and  Gautier,  C.:  Mesoscale  events  within  a  GATE  tropical  depression.  Mon. 
Wea.  Rev.,  106,  789-805,  1978. 

1255 


1260 


1265 


55 


Exp. 

no 

radiation 

radiation 

moist 

dry 

small 

large 

weak 

strong 

SST29 

SST28 

SST27 

2 

© 

© 

© 

© 

© 

3 

© 

© 

© 

© 

© 

4 

© 

© 

© 

© 

© 

5 

© 

© 

© 

© 

© 

6 

© 

© 

© 

© 

© 

7 

© 

© 

© 

© 

© 

8 

© 

© 

© 

© 

© 

9 

© 

© 

© 

© 

© 

10 

© 

© 

© 

© 

© 

11 

© 

© 

© 

© 

© 

12 

© 

© 

© 

© 

© 

13 

© 

© 

© 

© 

© 

14 

© 

© 

© 

© 

© 

15 

© 

© 

© 

© 

© 

16 

© 

© 

© 

© 

© 

1270 

Table  1,  Experiments  2-16  categorized  by  either:  no  radiation,  or  radiation  included; 
moist,  or  dry;  small,  or  large;  weak,  or  strong;  value  of  SST,  either:  29,  28,  or  27  °C. 
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TABLE  2.  General  statistics  for  the  experiments:  Shown  are  the  time  the  maximum 
1305  azimuthally-averaged  tangential  wind  speeds  near  the  surface  reach  12  m  s'1;  the  radius 
of  maximum  winds  (RMW)  at  this  time;  the  pathway  taken  to  genesis;  the  maximum 
surface  tangential  winds  prior  to  the  SSCV  forming  for  systems  that  develop  along 
pathway  two;  the  time  at  which  the  system  becomes  a  tropical  storm;  the  RMW  at  this 
time;  the  time  the  system  becomes  a  hurricane;  the  RWM  at  this  time. 
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TABLE  3.  Properties  of  the  prominent  mid-level  vortex  (MLV)  that  forms  for  systems 
that  develop  along  pathway  two:  Shown  are  the  maximum  azimuthally  averaged 
tangential  velocities  of  the  MLV,  its  RMW  and  height  at  the  time  the  SSCV  forms;  the 
1320  time  the  SSCV  forms;  the  maximum  surface  tangential  winds  and  the  RMW  at  the 
surface,  at  this  time. 
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Figure  Captions 

1340 

Figure  1.  Initial  conditions  for  a  small  weak  vortex  with  a  moisture  anomaly,  (a) 
Tangential  velocity.  The  contour  interval  is  1  ms'1,  (b)  Potential  temperature  anomaly. 
The  contour  interval  is  0.25  K.  (c)  Water  vapor  mixing  ratio.  The  contour  interval  is  1.5  g 
kg1. 

1345 

Figure  2.  Time  series  for  Experiment  1,  the  no  ice  case,  (a)  The  minimum  surface 
pressure,  (b)  The  azimuthally-averaged  maximum  tangential  velocity  and  the  height  it 
occurs,  (c)  The  azimuthally-averaged  tangential  velocity  at  the  lowest  model  level  and 
the  radius  it  occurs. 

1350 

Figure  3.  As  in  Figure  2,  for  Experiment  2. 

Figure  4.  Azimuthally-averaged  tangential  velocity  V(r,z)  for  Experiment  2,  at  (a)  t=79 
h,  and  (b)  t=83  h.  The  contour  interval  is  1  m  s'1. 

1355 

Figure  5.  As  in  Figure  2,  for  Experiment  3. 

Figure  6.  The  average  relative  humidity  for  a  circular  area  with  radius  of  100  km  from 
the  vortex  center,  at  levels  z=2.6,  4.9  and  7.3  km,  for  Experiment  3. 

1360 

Figure  7.  As  in  Figure  2,  for  Experiment  5. 

Figure  8.  As  in  Figure  2,  for  Experiment  6. 

1365  Figure  9.  As  in  Figure  2,  for  Experiment  7. 

Figure  10.  Time  series  for  Experiment  15.  (a)  Minimum  surface  pressure,  including 
comparison  with  the  non-radiation  case,  (b)  The  azimuthally-averaged  maximum 
tangential  velocity  and  the  height  it  occurs,  (c)  The  azimuthally-averaged  tangential 
1370  velocity  at  the  lowest  model  level  and  the  radius  it  occurs,  (d)  Total  liquid  and  ice  mass 
content  in  a  cylindrical  volume  of  radius  100  km  from  the  center,  (e)  Clear  sky  short 
wave  radiation  at  the  surface. 

Figure.  11.  Azimuthally-averaged  fields  of  (a)  Mixing  ratio  of  ice  (g  kg'1),  and  (b)  radial 
1375  component  of  wind  velocity  (m  s'1),  at  t=75  h,  for  Experiment  2. 

Figure.  12.  Comparison  of  Experiments  2  and  5.  (a)  The  total  mass  of  ice  in  the  fine  grid 
domain,  (b)  The  total  mass  of  ice  within  a  radius  of  100  km  of  the  center  of  the  domain, 
(d)  The  vertical  mass  transport  of  air  within  a  radius  of  100  km  of  the  center  of  the 
1380  domain,  at  z=7.35  km. 

Figure.  13.  The  azimuthally-averaged  potential  temperature  perturbation  for  Experiment 
2,  at  t=79  h.  The  contour  interval  is  0.5  K. 
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1385  Figure.  14.  Horizontal  sections  for  Experiment  2,  at  t=78  h.  (a)  Vertical  vorticity  (xlO’4 
rad  s’1),  (b)  potential  temperature  (K),  (c)  pressure  (mb),  at  z=98  m,  and  (d)  potential 
temperature  (K),  at  z=7.35  km. 

Figure.  15.  Horizontal  sections  for  Experiment  2,  at  t=80  h.  (a)  Vertical  vorticity  (xlO’4 
1390  rad  s’1),  (b)  potential  temperature  (K),  (c)  water  vapor  mixing  ratio  (g  kg'1),  (d)  pressure 
(mb),  at  z=98  m,  and  (d)  vertical  velocity  (m  s’1),  at  z=1.6  km. 

Figure.  16.  Time  series  of  the  average  potential  temperature  and  water  vapor  mixing 
ratio  within  100  km  of  the  center,  at  z=305  m,  for  Experiment  2. 
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Figure  1.  Initial  conditions  for  a  small  weak  vortex  with  a  moisture  anomaly,  (a)  Tangential 
velocity.  The  contour  interval  is  1  m  s'1,  (b)  Potential  temperature  anomaly.  The  contour 
interval  is  0.25  K.  (c)  Water  vapor  mixing  ratio.  The  contour  interval  is  1.5  g  kg'1. 
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Figure  2.  Time  series  for  Experiment  1,  the  no  ice  case,  (a)  The  minimum  surface 
pressure,  (b)  The  azimuthally-averaged  maximum  tangential  velocity  and  the  height  it 
1440  occurs,  (c)  The  azimuthally-averaged  tangential  velocity  at  the  lowest  model  level  and  the 
radius  it  occurs. 
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Figure  3.  As  in  Figure  2,  for  Experiment  2. 
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Figure  4.  Azimuthally-averaged  tangential  velocity  V(r,z)  for  Experiment  2, 
at  (a)  t=79  h,  and  (b)  t=83  h.  The  contour  interval  is  1  m  s'1. 
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Figure  5.  As  in  Figure  2,  for  Experiment  3. 
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Figure  6.  The  average  relative  humidity  for  a  circular  area 
with  radius  of  100  km  from  the  vortex  center,  at  levels 
z=2.6,  4.9  and  7.3  km,  for  Experiment  3. 
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Figure  7.  As  in  Figure  2,  for  Experiment  5. 
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Figure  8.  As  in  Figure  2,  for  Experiment  6. 
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Figure  9.  As  in  Figure  2,  for  Experiment  7. 
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Figure  10.  Time  series  for  Experiment  15.  (a)  Minimum  surface  pressure,  including 
comparison  with  the  non-radiation  case,  (b)  The  azimuthally-averaged  maximum 
tangential  velocity  and  the  height  it  occurs,  (c)  The  azimuthally-averaged  tangential 
velocity  at  the  lowest  model  level  and  the  radius  it  occurs,  (d)  Total  liquid  and  ice 
mass  content  in  a  cylindrical  volume  of  radius  100  km  from  the  center,  (e)  Clear  sky 
short  wave  radiation  at  the  surface. 
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Figure.  11.  Azimuthally-averaged  fields  of  (a)  Mixing  ratio  of  ice  (g  kg"1),  and 
(b)  radial  component  of  wind  velocity  (m  s"1),  at  t=75  h,  for  Experiment  2. 
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Figure.  12.  Comparison  of  Experiments  2  and  5.  (a)  The  total  mass  of  ice  in 
the  fine  grid  domain,  (b)  The  total  mass  of  ice  within  a  radius  of  100  km  of 
the  center  of  the  domain,  (d)  The  vertical  mass  transport  of  air  within  a 
radius  of  100  km  of  the  center  of  the  domain,  at  z=7.35  km. 
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Figure.  13.  The  azimuthally-averaged  potential  temperature 
perturbation  for  Experiment  2,  at  t=79  h.  The  contour  interval  is 


0.5  K. 
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Figure.  14.  Florizontal  sections  for  Experiment  2,  at  t=78  h.  (a)  Vertical  vorticity  (xl0“ 
rad  s’1),  (b)  potential  temperature  (K),  (c)  pressure  (mb),  at  z=98  m,  and  (d)  potential 
temperature  (K),  at  z=7.35  km. 
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igure.  15.  Horizontal  sections  for  Experiment  2,  at  t=80  h.  (a)  Vertical  vorticity  (xl(T 
potential  temperature  (K),  (c)  water  vapor  mixing  ratio  (g  kg'1),  (d)  pressure  (mb),  at  z= 
(d)  vertical  velocity  (m  s'1),  at  z=1.6  km. 
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Figure.  16.  Time  series  of  the  average  potential  temperature  and  water  vapor  mixing 
ratio  within  100  km  of  the  center,  at  z=305  m,  for  Experiment  2. 
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